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ABSTRACT 

(Distribution  Limitation  Statement  No.  2) 


This  report  presents  the  results  of  an  analytical  research  effort  aimed  toward 
defining  the  degradation  of  structural  integrity  that  may  result  when  design 
conditions  are  not  met  during  th°  construction  of  hardened  facilities.  Param¬ 
eters  of  primary  concern  were:  (a)  variations  in  the  soil  resistance  around  a 
buried  structure  (nonuniform  backfill  properties);  (b)  variations  in  the  struc¬ 
tural  geometr,  (improperly  enforced  construction  tolerances);  (c)  variations  in 
structural  material  properties.  A  computer  program  was  developed  to  study  the 
influence  of  variations  in  these  parameters  by  representing  the  structure-soil 
system  by  a  lumped-parameter  model  and  solving  the  resulting  equations  by  numer¬ 
ical  integration.  The  program  is  restricted  to  planar  reinforced  concrete 
structures,  with  particular  emphasis  on  buried  arches  and  cylindrical  structures. 
The  program  considers  both  elastic  and  inelastic  response  of  the  structure  under 
anv  combination  of  axial  force  and  bending.  It  is  also  capable  of  representing 
a  variety  of  multilinear  blast  pressure  pulses  that  may  approach  the  structure 
from  any  direction  in  the  plane  of  the  structure. 
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SECTION  I 


INTRODUCTION 


1 .  Objectives: 

The  basis  of  the  problem  considered  in  this  study  is  an  acknowledgment 
of  the  fact  that  significant  differences  may  exist  between  the  conditions 
assumed  for  the  design  of  hardened  underground  structures  and  the  correspond¬ 
ing  conditions  that  actually  exist  in  the  as-built  structure.  Of  particular 
concern  are  the  variations  that  might  exist  between  the  assumed  and  the  real 
soil  resistance  characterist ics ,  the  geometry  or  alignment  of  the  structure, 
and  the  properties  of  the  structural  materials.  Hence,  the  primary  objective 
of  this  study  was  to  develop  a  method  to  investigate  the  influence  of  variations 
in  these  several  significant  parameters  on  the  resistance  of  a  burled  structure 
to  nuclear  blast  loading. 

A  secondary  objective  of  the  project  was  to  use  the  method  of  analysis 
thus  developed  to  investigate  to  a  limited  extent  the  influence  of  variations 
in  some  of  the  more  significant  parameters  on  the  blast  resistance  of  buried 
structures. 


2 .  Scope  of  Effort: 

Because  a  comprehensive  study  of  the  problems  indicated  in  the  subject 
of  this  report  is  formidable,  the  scope  of  the  study  must  be  clearly  de¬ 
fined.  The  major  effort  was  to  develop  an  analytical  method  to  study  the 
influence  of  significant  parameter  variations  on  the  resistance  of  buried 
structures  to  nuclear  blast  loads.  A  secondary  effort  was  to  use  this  method 
to  study  the  influence  of  such  variations  on  specific  structures. 

It  is  emphasized  that  the  significance  of  parameter  variations  on  the 
resistance  of  buried  structures  is  the  primary  interest  of  this  study.  The 
research  reported  here  is  not  intended  as  a  further  study  of  soil-structure 
interaction  or  to  identify  quantitatively  the  influence  of  structural  deforma¬ 
tions  on  the  pressures  that  exist  across  a  soil-structure  interface. 

Because  of  its  wide  use  in  protective  construction  and  its  probable 
sensitivity  to  the  effects  of  parameter  variations,  the  buried  reinforced 
concrete  arch  was  used  as  the  basis  structural  type  in  this  study.  For  example, 
consider  a  fully  buried  arch  that  has  been  designed  on  the  assumption  that  the 
soil  resistance  around  it  is  uniform  at  a  given  level  of  resistance.  These 
studies  relate  to  the  extent  to  which  the  resistance  of  such  an  arch  is  reduced 
if  over  part  of  its  surface  area  the  soil  resistance  is  reduced  by  ten  percent, 
fifteen  percent,  thirty  percent,  etc,,  of  the  basic  resistance  offered  by  the 
soil  over  the  remaining  surface  area  c*  the  arch.  Similarly,  If  an  arch  Is 
designed  on  the  assumption  that  it  is  an  arc  of  a  true  circle,  by  what  percent¬ 
age  is  its  resistance  reduced  if  during  the  construction  and  backfilling  oper¬ 
ations  variations  from  the  assumed  circular  shape  equal  to  one  percent,  two 
percent,  five  percent, etc. , of  the  radius  are  developed? 


The  method  of  analysis  that  was  developed  is  sufficiently  general  that  it 
can  be  used  directly  for  the  study  of  any  two-dimensional  reinforced  concrete 
configuration  having  a  constant  cross  section.  Thus,  it  is  directly  applicable 
to  the  analysis  of  cylindrical  or  rectangular  reinforced  concrete  structures 
of  constant  wall  thickness.  With  relatively  minor  modifications,  the  method 
can  also  be  applied  to  other  planar  structures  of  other  materials. 

The  method  of  analysis  that  was  developed  was  a  computer  program  capable 
of  solving  for  the  response  of  a  lumped  parameter  model  of  the  soil-structure 
systems  to  blast- induced  soil  pressures.  Response  was  not  restricted  to  the 
linearly  elastic  range;  inelastic  structural  behavior  under  any  combination  of 
axial  force  and  bending  can  also  be  considered. 


3 .  Significance  of  Study : 

If  properly  exploited  the  computer  program  developed  in  this  study  could 
yield  results  of  substantial  immediate  significance  not  only  to  the  design  and 
construction  of  hardened  facilities  in  the  future  but  also  to  the  evaluation 
of  the  hardness  of  existing  similar  facilities.  For  example,  in  the  design  of 
future  facilities,  e  better  understanding  of  the  influence  of  variations  in 
parameters  such  as  soil  resistance,  structural  material  properties,  and  structural 
geometry  would  be  of  great  value  to  the  designer  as  he  prooortions  a  structure 
to  provide  a  specified  level  of  blast  resistance.  Considering  the  uncertainties 
inherent  In  these  parameters  in  any  design  problem,  a  proper  approach  is  to 
design  the  structure  to  provide  a  specified  probability  that  the  facility  will 
survive,  with  acceptable  damage  levels,  a  given  intensity  of  blast  loading. 

The  computer  program  developed  herein  can  be  used  effectively  to  obtain  statis¬ 
tical  data  on  the  influences  of  the  several  parameters  on  the  resistance  that 
would  be  required  for  any  useful  probability  analysis. 

Similarly,  it  could  be  used  to  provide  the  input  information  that  would 
be  required  for  a  probabilistic  analysis  of  the  blast  resistance  of  in-place 
structures.  It  should  be  possible,  for  example,  to  determine  for  a  given  in- 
place  facility  the  extent  to  which  the  structural  geometry,  the  structural 
material  properties,  and  the  soil  properties  differed  from  the  values  of  these 
parameters  that  were  used  in  the  design  of  the  facility.  On  the  basis  of  this 
information,  the  computer  program  could  be  used  to  determine  the  extent  to  which 
the  actual  blast  resistance  of  the  facility  differ*-  from  that  for  which  it  was 
originally  designed. 

Further,  data  that  can  be  acquired  through  the  exploitation  of  this  computer 
program  should  be  especially  useful  in  the  preparation  of  specifications  for  the 
construction  of  hardened  facilities  In  the  future.  A  better  understanding  of 
the  influences  of  variations  in  the  several  significant  parameters  on  the 
resistance  of  the  structure  is  essential  in  order  to  specify  permissible  design 
ar.d  construction  tolerances. 


SECTION  II 


METHOD  OF  STUDY 


1 •  Seneral  Discussion  of  Problem: 


a  .  ,An  actual  structure  will  inevitably  differ  from  the  planned  structure 
and  from  the  idealized  structure  that  was  considered  in  the  design.  Errors 
afinn*  c  1  rc^s tances  in  tie  construction  can  result  in  marked  varl- 

b^nMuf-  !  C°^  sPeclfled  the  Plans  which,  in  turn,  may  have 

been  s impl  if  ted  and  idealized  to  an  appreciable  extent  to  arrive  at  the 
design.  Variations  of  this  type  always  exist  in  any  structure  but  they 
take  on  higher  than  usual  importance  when  considered  in  regard  to  underground 
structures  that  are  designed  for  high  levels  of  nuclear  blast  loads.  9 


The  principal  differences  that  may  exist  between  the 
assumed  for  design  and  the  real  structure  as  built  In  its 
ment  can  be  grouped  as  follows: 


idealized  structure 
underground  envlron- 


(1)  Alignment  --  the  differences  between  the  actual  and  the  assumed 
locations  of  the  axes,  planes,  or  other  surfaces  that  define  the 
geometry  of  the  structural  member  or  element. 

(2)  Proportions  —  the  differences  between  the  ideal  and  the  actual 
dimensions  of  structural  elements;  for  example,  variations  In  the 
depth  of  a  reinforced  concrete  section,  the  amount,  spacing,  and 

position  of  the  reinforcement  in  the  slab,  and  similar  other  varla- 
t Ions . 


(3)  Structural  material  properties  -  differences  between  the  assumed 
wereusedCtUa  phySlcaI  ProPertIes  the  structural  materials  that 


(4)  Soil  resistance  properties  —  differences  between  the  ideal  and  the 
actual  physical  properties  of  the  soils  or  other  materials  that 
support  the  structure  and  within  which  the  structure  is  placed. 

This  includes  not  only  variations  in  the  basic  level  of  soil  re¬ 
sistance  assumed,  but  also  variations  in  the  magnitude  of  this 
resistance  as  it  may  differ  from  point  to  point  around  the  structure. 


The  variations  to  be  expected  in  the  several  primary  parameters  Identi¬ 
fied  above  derive  from  a  variety  of  sources  and  are  subject  to  varying  de- 

beeDo«ib?0rT01  ^  -hC  st!*ucfural  designer  and  the  builder.  Although  it  may 
e  possible  to  exercise  only  little  control  over  the  variations  that  will 
.occur  natural ly  in  many  structural  and  soil  material  properties,  it  is 

mannlJurfp  th*  influence  of  thcse  variations,  of  whatever 

magnitude  and  font,  they  may  take,  on  the  behavior  and  hence  on  the  ultimate 
resistance  of  the  structure  being  designed. 


Of  perhaps  even  greater  importance  is  the  need  for  a  better  understanding 
Sl!OUld  bC  undertaken  fn  the  construction  process  to  con- 
ha  tH  ?r,a- '0nS  th3t  Can  deveI°P  In  those  parameters  over  which  we  can 
have  an  effective  measure  of  control.  For  example,  it  is  possible  to  exercise 
g  er  levels  of  control  and  to  Insist  on  construction  to  closer  tolerances 
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(though  admittedly  with  an  increased  cost  factor),  if  it  can  be  demonstrated 
that  such  more  rigid  construction  specifications  are  essential  to  the  attain¬ 
ment  of  a  specified  level  of  protection.  Conversely,  if  It  cart  be  shown  that 
somewhat  greater  construction  tolerances  are  permissible  without  a  signifi¬ 
cant  reduction  in  the  resistance  level  of  the  resulting  structure,  then  the 
economies  associated  with  less  restrictive  specifications  can  be  achieved. 

In  any  event,  without  information  beyond  that  which  we  now  possess  in  regard 
to  the  significance  of  variations  in  these  several  parameters,  it  is  im¬ 
possible  to  prepare  a  set  of  construction  specifications  that  are  consistent 
with  the  physical  conditions  that  prevail  and  with  a  specified  probability  of 
achieving  a  stipulated  level  of  protection  to  a  given  set  of  input  loadinq 
conditions.  3 

All  the  parameters  identified  above  are  significant  in  the  full  evalua¬ 
tion  of  the  resistance  of  a  given  soil-structure  system;  but  the  influence  of 
variations  in  some  of  these  parameters  is  more  easily  dealt  with  than  are 
variations  in  others.  For  many  comparatively  simple  geometric  forms  and  loading 
conditions,  the  influence  of  variations  in  structural  material  properties  can 
be  easily  determined.  It  may  be  true  that  to  establish  the  significance  of 
these  variations  a  large  amount  of  more  or  less  routine  computational  effort 
will  be  required;  but  simple  methods  exist  for  the  study  of  the  influence  of 
such  parameter  variations.  Similarly,  methods  are  also  readily  available 
whereby  for  simple  geometrical  forms  the  influence  of  unintended  variations  in 
t  e  proportions  (widths,  thicknesses,  reinforcement  steel  percentages,  etc.) 
on  the  resistance  of  the  structural  element  are  also  easily  evaluated;  although 
again,  a  substantial  amount  of  routine  computation  effort  may  be  required. 

At  the  present  time,  however,  no  procedures  exist  whereby  the  effects  of 
departures  from  the  specified  structural  geometry  (variations  in  ladius 
around  an  arch,  for  example)  or  In  the  uniformity  of  soil  resistance  around 
the  structure  can  be  readily  determined.  Consequently,  in  this  effort  primary 
attention  was  given  to  the  study  of  the  influence  of  variations  In  these  latter 
two  parameters,  although  the  program  that  was  developed  is  capable  also  of 
considering  the  influence  of  variations  in  the  former  two  parameters  as  well. 


2 .  General  Method  of  Analysis: 

To  obtain  a  program  of  sufficiently  general  character  that  it  would  have 
wide  use  In  the  study  of  the  problem  described  above,  It  was  necessary  to  re¬ 
present  the  buried  structure  and  the  soil  which  surrounds  I t  as  a  lumped 
parameter,  soil-structure  model,  the  response  of  which  to  superimposed  blast 
loadings  could  be  studied  by  existing  numerical  methods  of  analysis.  Conse¬ 
quently,  the  soil-structure  system  was  represented  as  a  series  of  concentrated 
masses  which  were  connected  by  massless  bars;  these  bars,  though  massless,  re¬ 
tained  their  axial  and  bending  stiffnesses.  Similarly,  the  applied  blast- 
induced  loads,  the  soil  resisting  forces,  and  the  internally  generated  struc¬ 
tural  resistances  were  also  concentrated  at  the  several  mass  points. 

Although  they  will  be  discussed  in  detail  in  later  sections  of  this 
report,  it  is  appropriate  to  identify  here  the  elements  of  the  mathematical 
model  to  which  attention  must  be  given.  The  problem  areas  encountered  in  the 
development  of  a  lumped  parameter  model  that  would  be  equivalent  to  the  real 
system  include  the  following: 
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(0  The  number  and  magnitude  of  the  masses  used  to  represent  the  system. 

(2)  Evaluation  of  the  internal  structure  resistances,  and  the  concen- 
t rat  ion  of  these  resistances  at  the  mass  points, 

(3)  Representation,  as  a  function  of  displacement,  of  the  resistance 
offered  by  the  soil  to  motion  of  the  structure,  and  the  concentra¬ 
tion  of  these  soil  resisting  forces  at  the  mass  points. 

(4)  Representation  of  the  blast- induced  free-fleld  forces  in  the  soil 
surrounding  the  structure,  and  the  concentration  of  these  blast 
forces  at  the  mass  points  of  the  model. 

(5)  Establishment  of  criteria  for  ''failure"  of  the  structure  to  serve 

as  a  basis  for  comparison  of  the  effects  of  variations  In  the  several 
parameters  of  interest. 


.  .  h?Wn9  defined  the  several  elements  of  the  mathematical  model 

identified  above,  it  was  then  possible  to  write  a  computer  program,  using 
numerical  methods,  to  solve  for  the  response  of  the  model  to  a  variety  of  load 
functions  that  may  be  imposed  upon  it.  A  detailed  description  of  the  com¬ 
puter  program  that  was  developed  for  this  purpose  is  given  in  Section  III  of 
tn is  report. 


3 •  Development  of  Lumped  Mass  Model; 

eut.  1,16  sev®*'al  of  the  mathematical  model  of  the  actual  sol  1-structun 

system  were  identified  in  the  preceding  section.  In  the  paragraphs  that  follow 

g°vendeta  1Cd  dlscuss  on  of  the  Pr°Per  evaluation  of  each  of  these  elements  Is 

l  -^er  M*ss  po»nts:  It  is  impossible  to  specify  in  general  terms 
the  number  of  mass  points  that  should  be  selected  to  represent  a  particular 
real  system.  The  number  of  mass  points  that  should  be  used  in  a  given  case 
depends  upon  the  geometrical  complexity  of  the  system  being  represented,  the 
egree  of  symmetry  in  the  system,  and  the  accuracy  to  be  sought  in  the  solution. 
Consequently,  the  computer  prodram  was  prepared  to  permit  the  use  of  any  number 
o  mass  points.  As  far  as  the  program  itself  is  concerned,  the  number  of  mass 
points  used  is  limited  only  by  the  size  of  the  available  computer  system. 

.  illS-Pl  Equivalent  Masses:  For  a  structural  element  for  which  the 

dynamic  response  to  a  given  loading  is  known  or  can  be  determined,  it  is 
possible  to  define  the  size  of  the  masses  and  the  character  of  the  resistance 
functions  that  are  required  to  ensure  equivalence  between  the  response  of  the 
structural  element  and  the  response  of  a  lumped  parameter  model  of  it.  It  Is, 
owever,  impossible  to  do  this  for  the  very  complex  soil-structure  system  that 
s  being  considered  here.  There  is  some  uncertainty  In  regard  to  the  portion 
o  the  structure  that  should  be  lumped  at  the  mass  concentration  points;  but 
“ncertainty  exists  in  regard  to  the  amount  of  the  surrounding  soil 
that  should  be  assumed  to  respond  with  the  structure  and,  hence,  that  should 
be  included  in  the  concentrate  i  mass  of  the  mathematical  model. 
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1  2 

Some  guidance  In  this  regard  Is  given  by  Whipple  ’  from  whose  work  It 
may  be  concluded  that  the  period  of  vibration  of  a  burled  arch,  as  it  responds 
to  the  forces  imposed  by  the  passage  of  an  airblast  wave  across  the  ground 
surface  above  it,  is  probably  about  the  same  as  the  period  of  the  same  arch 
when  completely  above  ground.  That  is  to  say,  the  mass  of  the  responding  system 
should  be  increased  because  of  the  surrounding  soil  that  moves  with  it  by 
approximately  the  same  amount  by  which  the  resistance  of  the  arch  is  augmented 
by  che  resistance  of  the  soil  adjacent  to  it.  Consequently,  assuming  the 
structure  to  be  such  that  its  periods  of  vibration  can  be  evaluated  in  an 
aboveground  situation,  the  equivalent  mass  and  soil  resistance  that  should  be 
employed  In  the  represent?  i.  Ion  of  the  buried  soil-structure  system  could  be 
determined  by  trial  and  error  to  give  approximately  the  same  periods  of  vibra¬ 
tion  that  would  exist  in  the  aboveground  case. 

Some  further  insight  into  this  question  is  given  by  the  studies  which 
show  that  for  a  perfect  fluid,  a  volume  of  fluid  included  within  a  distance  of 
one  radius  of  a  submerged  cylindrical  section  can  be  assumed  to  respond  with 
the  circular  structure  and  hence  should  be  considered  as  effective  mass.  Soil, 
however,  Is  obviously  not  a  fluid;  hence,  the  extent  to  which  this  criterion 
for  effective  mass  is  applicable  is  uncertain. 

In  view  of  these  factors,  the  computer  program  for  the  analysis  of  the 
structures  to  be  considered  here  permits  the  use  of  any  effective  mass;  it 
remains  for  the  user  of  the  program  to  select  the  effective  mass  that  he  con¬ 
siders  most  appropriate  for  the  conditions  being  studied.  It  should  also  be 
noted  that  for  the  parameter  studies  that  will  be  made  using  this  program,  tie 
question  of  the  correct  effective  mass  is  probably  not  important  in  a  given 
case.  It  will  be  remembered  that  the  objective  of  the  project  was  to  develop 
a  program  to  study  the  influence  of  variations  in  parameters;  the  magnitude 
of  the  soil  mass  can  itself  be  considered  as  such  a  parameter  and  the  signifi¬ 
cance  of  variations  in  it  can  be  evaluated  in  this  fashion.  If  the  results  of 
such  studies  should  indicate  that  the  amount  of  soil  that  is  assumed  to  re¬ 
spond  with  the  buried  structure  influences  substantially  the  structure's 
resistance  to  blast  forces,  this  is  in  itself  a  significant  element  of  infor¬ 
mation  and  one  that  would  indicate  strongly  the  need  for  a  research  effort 
directed  specifically  to  this  problem. 

(c)  Representation  of  Soil  Resistance:  Because  of  the  nature  of  this 
study,  no  attempt  was  made  to  relate  the  soil  resistance  function  used  in  the 
mathematical  model  to  any  real  soil  situation.  It  was,  however,  necessary 
that  the  basic  soil  resistances  be  reasonable  and  that  provisions  be  made  in 
the  computer  program  not  only  to  study  the  influence  on  the  structural  re¬ 
sistance  of  variations  in  the  basic  resistance,  but  also  to  permit  the  study 
of  the  influence  of  a  variable  soil  resistance  around  the  structure.  It  will 


^Whipple,  C.  R.,  "The  Dynamic  Response  of  Shallow-Buried  Arches  Subjected  to 
Blast  Loading,"  Report  for  Air  Force  Special  Weapons  Center,  Contract  No. 

AF  29(601)-2591,  Project  No.  1080,  University  of  Illinois,  1961.  (UNCLASSIFIED) . 

2 

'Whipple,  C.  R.,  "Numerical  Studies  of  the  Dynamic  Response  of  Shallow  Buried 
Arches  Subjected  to  Blast  Loading,"  Draft  of  Report  for  Air  Force  Special 
Weapons  C.encer,  Contract  No.  AF  29 (601 ) -259 1  and  AF  29(601)-4508,  Project  No. 
1080;  University  of  Illinois,  1961  (UNCLASSIFIED). 
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be  remembered  that  one  of  the  parameter  to  which  attention  was  to  be  direct‘d 
initially  was  the  Influence  of  the  introduction  of  ‘’soft  spots"  at  various 
points  around  the  structure  in  relation  to  the  level  of  soil  resistance  that 
exists  generally  around  the  structure. 

It  Is  recognized  t.iat  in  the  real  case  the  soil  resistance  to  motion  of 
a  structure  embedded  within  it  is  a  highly  variable  and  uncertain  function, 
and  that  it  depends  on  the  soil  type,  the  moisture  conditions,  the  in  situ 
stress  conditions,  the  nature  and  distribution  of  the  motions  of  the  deforming 
structure,  and  many  other  factors.  It  was  dec!dru,  however,  that  for  this 
initial  study  the  soil  resistance  would  be  restricted  to  the  case  of  a  linear! 
elastic  soil  which  was  capable  of  carrying  only  compressive  stress,  and  which 
unloaded  to  zero  stress  on  the  same  modulus  of  deformation  upon  which  it  was 
loaded. 

Although  this  form  of  resistance  function  is  obviously  not  a  true  re¬ 
presentation  of  any  actual  soil,  it  is  considered  to  be  a  reasonable  approx¬ 
imation  of  the  behavior  of  many  soil  types,  at  least  under  conditions  of  re¬ 
latively  small  deformation.  Before  more  extensive  studies  are  made  into  the 
significance  of  a  wide  variety  of  soil  resistance  forms,  it  seems  appropriate 
that  a  more  nearly  complete  understanding  be  acquired  of  the  effects  of  vari¬ 
ations  in  soil  resistances  of  this  relatively  simple  and  basic  form.  The 
computer  program  has,  therefore,  been  written  for  this  type  of  soil  resistance 
function;  but  it  could  be  modified  without  great  effort  to  permit  the  study 
of  more  complex  load  forms.  A  further  discussion  of  the  soil  resistance  func¬ 
tions  and  the  way  in  which  the  soil  resistance  was  concentrated  at  the  several 
mass  points  Is  given  In  the  detailed  discussion  of  the  development  of  the 
computer  program  in  Section  III* 

(d)  Representation  of  Load  Function:  From  the  standpoint  of  representing 
the  actual  soil-structure  system  by  an  equivalent  lumped-mass  mathematical 
model,  the  proper  representation  of  the  blast- I nduced  load  function  is  the  most 
difficult  problem  and  is  probably  the  source  of  greatest  uncertainty.  In  the 
actual  physical  case,  it  is  clearly  impossible  to  separate  the  soil  resistance 
function  and  the  applied  load  funtion  on  the  structure  for  the  obvious  reason 
that  these  two  factors  combine  in  some  fashion  to  prod  ice  the  pressure  that 
exists  on  the  soil-structure  interface.  This  question  is,  of  course-  the  very 
heart  of  the  soil-structure  interaction  problem  which  is  receiving  much  re¬ 
search  emphasis,  but  which,  as  stated  at  the  outset,  is  beyond  the  scope  of 
this  present  study. 

It  is  clearly  incorrect  in  the  real  case  to  assume  that  the  pressures 
transmitted  to  the  buried  structure  from  the  b! ast- induced  free-ficld  stresses 
in  the  soil  surrounding  the  structure  are  independent  of  and  are  unaffected 
by  the  deformation  of  the  structure.  An  effort  was  made  in  the  early  stages 
of  this  investigation  to  recognize  the  interdependence  of  soil  resistance  and 
load  application  to  the  structure  by  endeavoriny  to  introduce  the  load  func¬ 
tions  as  displacements  of  the  bases  of  the  springs  that  were  used  in  the  model 
to  represent  the  soil  resistance.  Because  of  the  relationships  that  are 
known  to  exist  between  free-field  soil  stress  ant  particle  velocity  in  the  soil, 
this  approach  to  the  soil-structure  interaction  problem  is  still  thought  to 
be  potentially  valid.  However,  attempts  to  arrive  at  a  satisfactory  model 
using  this  approach  have  not  as  yet  been  successful;  the  primary  source  of 
difficulty  is  an  inability  to  represent  properly  the  load  function  as  it  is 
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influenced  by  the  relative  deformations  of  the  soil  spring  base  and  of  the 
structure  at  its  point  of  attachment  to  the  spring.  It  is  a  relatively 
simple  matter  to  translate  a  free-field  soil  stress  pulse  into  corresponding 
free-field  solid  particle  motions  and  hence  into  resulting  forces  in  a  soil 
spring.  To  account  at  the  same  time  for  the  effect  on  the  spring  force  of 
the  motion  of  the  structure,  which  is  a  combination  of  the  motion  of  the 
structure  relative  to  the  soil  immediately  around  it  as  well  as  of  a  rigid 
body  motion  of  the  structure  as  it  moves  in  some  manner  with  the  soil  surround¬ 
ing  it,  is  much  more  difficult. 

Therefore,  acknowledging  the  incorrectness  of  it,  but  recognizing  at  the 
same  time  that  there  presently  exists  no  better  alternative,  we  assumed  the 
blast- induced  load  functions  used  for  this  study  to  be  directly  related 
to  and  derivable  from  the  free-field  blast- induced  stresses  that  exist  in 
the  soil  surrounding  the  structure.  It  is  appropriate  to  note  here  that  the 
computer  program  was  written  so  that  a  variety  of  multilinear  free-field  soil 
pressure  functions  could  be  studied;  the  details  of  these  pressure  function 
representations,  and  the  techniques  that  were  employed  to  accumulate  the  free- 
field  pressures  and  concentrate  them  as  forces  applied  directly  to  the  mass 
points  of  the  mathematical  model  are  contained  in  Section  III. 

One  further  restriction  in  regard  to  the  load  pulse  should  also  be  noted. 
No  provision  has  been  made  for  varying  the  shape  of  l.'e  load  p  ‘>e  as  it 
propagates  across  the  structure;  the  same  basic  free-field  stress  pulse  is 
assumed  to  exist  at  all  points  in  the  soil  which  surrounds  the  structure. 
However,  the  program  is  written  so  that  the  pressure  pulse  may  approach  the 
structure  and  traverse  it  in  any  direction  that  may  be  appropriate. 

(e)  Structural  Resistance:  As  noted  earlier,  the  structure  is  repre¬ 
sented,  for  purposes  of  analysis,  as  a  series  of  concentrated  masses  connected 
by  massless  bars  or  links  which,  though  massless,  develop  internal  resisting 
forces  as  deformations  are  imposed  upon  them.  Thus,  it  was  necessary  to 
develop  a  method  whereby  these  internal  resisting  forces,  which  act  as  reactions 
on  the  concentrated  mass  points,  may  be  computed  from  the  deformations  that 
are  imposed  upon  the  bars. 

In  the  process  of  the  dynamic  analysis  of  the  entire  lumped-mass  system, 
the  changes  in  the  coordinates  of  the  mass  points,  and  hence  of  the  ends  of 
each  connecting  bar,  will  be  determined.  From  these  changes  in  the  mass  point 
coordinates,  the  changes  in  connecting  bar  lengths  and  the  rotations  of  one 
bar  relative  to  the  adjacent  bars  can  be  readily  computed.  From  the  total 
cnange  in  bar  length,  which  is  assumed  to  be  uniformly  distributed  throughout 
the  length  of  the  bar,  and  from  the  total  rotation  of  adjacent  bars,  which  is 
assumed  to  be  the  accumulation  of  uniform  curvature  over  a  length  equal  to  the 
length  of  one  bar,  the  strain  distributions  that  will  exist  on  the  structural 
cross  sections  at  each  end  of  a  bar,  and  hence  on  each  side  of  a  mass  point, 
are  readily  determined.  Similarly  then,  from  these  strain  distributions,  the 
resulting  thrusts,  moments,  and  corresponding  shears  are  computed.  The  shears 
and  thrusts  thus  determined  are  the  internal  resisting  forces  that  act  on  the 
concentrated  mass  points. 

In  concept,  the  procedure  outlined  above  is  relatively  simple  and  straight¬ 
forward;  in  reality,  the  development  of  a  computer  program  to  perform  this 
phase  of  the  analysis  is  quite  complex.  The  program  must  be  sufficiently 
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general  to  evaluate  the$e  forces  for  a  reinforced  concrete  section  as  it  re- 
sponds  to  the  combined  effects  of  axial  stress  and  bending  well  beyond  the 
region  of  elastic  behavior  of  these  materials.  The  details  of  the  proqram 

toT^Pute  these  internal  resisting  structural  forces  are 
presented  in  Section  III. 

(0  Failure  Criteria:  Since  the  objective  of  this  study  was  to  deter¬ 
mine  the  influences  on  the  failure  load  of  a  burled  structure  of  variations  in 
significant  parameters  such  as  soil  resistance,  structural  geometry,  and 
material  properties,  it  is  necessary  to  define  a  structural  failure  condition 

to  serve  as  a  basis  against  which  these  comparative  parane  ter  effects  can  be 
measured . 


In  an  actual  case,  structural  failure  would  generally  be  defined  In 
terms  of  structural  deformations  or  motions  that  are  of  sufficient  magnitude 
to  impair. to  a  specified  degree  the  ability  of  the  structure  to  function  as 
it  was  originally  intended.  Since  this  study  is  not  related  to  a  particular 
structural  type  or  to  structures  designed  to  perform  a  particular  function,  it 
is  impossible  to  define  structural  failure  on  the  bases  outlined  above.  By 
making  comparatively  minor  changes  in  the  computer  program,  It  would  be 
possible  to  specify  a  variety  of  failure  criteria  as  they  may  be  appropriate 
for  various  structural  type-  and  conditions.  However,  in  the  program  as 
currently  writter,  failure  was  somewhat  arbitrarily  defined  to  exist  when  the 
deformations  in  the  structure  had  developed  in  sufficient  magnitude  at  enough 
points  in  tne  structure  to  transform  the  structure  into  a  mechanism.  Thus, 
for  a  two-hingea  arch  to  which  the  program  has  thus  far  been  applied,  failure 
was  assumed  to  have  occurred  when  general  yielding  developed  across  the  struc¬ 
tural  cross  section  at  two  points  on  the  circumference  of  the  arch. 


Although  a  mechanism  might  not  necessarily  be  developed,  it  was  also 
assumed  that  failure  would  have  occurred  when  the  concrete  strains  at  any  two 
points  in  the  arch  exceeded  the  concrete  strain  at  which  the  concrete  would 
have  crushed  and,  hence,  would  have  suffered  a  reduction  in  stress  with  further 
increase  in  strain.  It  was  further  assumed  that  a  combination  of  the  general 
yielding  case  and  the  concrete  crushing  case  at  any  two  points  in  the  arch  would 
also  constitute  failure.  Although  such  a  failure  criterion  may  well  be  in¬ 
appropriate  for  a  given  physical  situation,  it  is  considered  to  be  reasonable 
as  a  basis  upon  which  to  make  the  comparative  parameter  variation  effect  studies 
ror  which  this  computer  program  was  written. 
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SECTION  III 


DEVELOPMENT  OF  COMPUTER  PROGRAM 


1 •  General  Description  of  Program : 

The  program  was  developed  for  use  on  the  University  of  Illinois  IBM  7094 
computer.  It  was  written  in  FORTRAN,  but  it  was  compiled  and  run  on  FASTRAN, 
a  compiler  developed  by  the  Department  of  Computer  Science  at  the  University. 
FASTRAN  is  faster  and  mo  •«*  efficient  than  FORTRAN;  however,  in  almost  all 
other  aspects  it  is  the  rune  as  IBM's  FORTRAN. 

The  program  itself  consists  of  a  main  program  which  merely  handles  data, 
stops  the  computer  at  the  p-oper  time,  and  calls  several  subroutines  as  they 
are  needed.  A  total  of  eight  subroutines  was  used,  although  not  all  of  then 
are  called  from  the  main  program;  some  are  called  from  other  subroutines.  The 
eight  subroutines,  BETA,  CHECK,  DEFORM,  ETC  I TE,  EXTERN,  INIT,  LOCATE,  and  YIELD, 
each  of  which  will  be  discussed  in  detail  later,  perform  all  computations  re¬ 
quired  in  the  analysis. 

During  the  development  of  the  program,  each  of  the  subroutines  was 
written  and  checked  for  errors  before  proceeding  to  the  next.  As  a  result 
most  of  the  variables  used  are  in  DIMENSION  and/or  COMMON  STATEMENTS. 

A  logic  flow  diagram  ©fthe  analysis  used  is  shown  in  Figure  1.  This 
method  of  analysis  is  typical  of  familar  numerical  integration  techniques  of 
solving  vibration  problems;  hence,  detailed  discussion  of  it  he:  e  is  un¬ 
necessary.  However,  further  comment  in  regard  to  the  more  complex  aspects  of 
it,  as  applied  to  the  type  of  problem  being  considered  In  this  study,  is  given 
i.i  the  descriptions  of  the  several  subroutines  that  follow  ii  later  sections. 

As  discussed  earlier,  the  soil-stfuctute  system  is  approximated  as  a 
lumped-mass  model  consisting  of  masses,  bars,  and  soil  springs.  I'uch  a  model 
is  illustrated  in  Figure  2. 

The  main  program  reads  all  input  data  and  writes  it  out  in  the  form  of 
an  echo  print  to  facilitate  checking  the  data.  These  input  data  and  the  form 
in  which  they  are  entered  into  the  computer  are  presented  in  the  appendix. 

These  data  include  the  initial  geometry  of  the  model,  including  the  soil  mass 
that  is  assumed  to  respond  with  the  model;  a  description  of  the  stress-strain 
curves  for  concrete  and  for  steel;  a  description  of  the  blast  pressure  pulse, 
Including  its  orientation  with  respect  to  the  horizontal,  and  the  ratio  be¬ 
tween  horizontal  and  vertical  free-field  blast  pressure;  the  seismic  velocity 
of  the  soil,  and  the  soil  resistance  function;  the  number  of  fibers  used  to 
represent  the  cross  section  of  the  structure,  the  area  of  reinforcing  steel, 
and  the  location  of  this  steel  In  the  cross  section;  the  initial  time  for  the 
program;  and  several  constants  that  define  acceptable  limits  of  computational 
accuracy  and  control  the  print-out  from  the  program. 

The  geometry  of  the  model  is  read  into  the  computer  by  specifying  the 
total  number  of  joints,  or  mass  points,  in  the  structure  (limited  in  the 
program  as  written  to  eleven,  but  readily  changeable  to  any  desired  number), 


10 


Hf»rt 
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JtIL 
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Locata  mv*  front  on  structure 
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As  tu—  accelerat  Ions 


Compute  velocity  and  coordinate  changes 
from  assumed  accelerations 


Compute  spring  def onset  Ions  froai 
coordinate  chenges _ 


X 


Compute  blast  and  spring 

forces  on  masses 


Use  computed  accelerations 

as  assumed  accelerat Ions 
for  this  trial 


increment  time  by  At 


Compute  bar  deformations  and  rotet Ions 
_ from  coordinate  cha nga < 


X 


Compute  bar  strains  from  defonaatlons 


X 


Compute  thrusts,  moments 

fttma  frw  ftrtl.ni 


X 


Sum  all  force  components 

on  masses 


X 


Acceleration  -  Force/Mass 


NO 


X 


Do  computed  accelerations  agree 
mltS  assumed  accelerations? 


YES 


Write  answers  -  final  accelerations, 
_ velocities,  forces,  etc. _ 


I 


^^-<(Has  structure  f al  iad?^>- 


J8L 


Wrlte  failure 

3= 


ENO 


FIGURE  1.  LOGIC  FLOW  DIAGRAM 


FIGURE  2.  LUMPED  MASS  MODEL  OF  STRUCTURE 


the  length  of  the  bars  connecting  the  mass  points  (initially,  all  bars  have 
the  same  length.),  and  the  angle  that  each  bar  makes  with  the  horizontal. 

The  program  sets  the  x  and  y  coordinates  of  Joint  1  equal  to  zero,  and  then, 
in  the  main  program,  computes  the  coordinates  of  the  other  joints  from  the 
following  equations. 


x(i+l)  *  xi  +  L  C0Sa\ 
v(l+l)  -  V|  +  1  !'"“i 

in  which  L  Is  the  original  length  of  the  bars  In  the  model  and  a,  Is  the  acute 
angle  that  bar  i  makes  initially  with  the  horizontal.  ' 

After  the  input  data  have  been  read  in  and  the  mass  point  coordinates 
have  been  computed  as  described  above,  the  subroutine  INIT  is  called  to  initial¬ 
ize  the  several  variables  at  their  proper  values  and  the  subroutine  LOCATE  is 
called  to  place  the  pressure  pulse  on  the  structure  and  begin  the  actual 
analysis. 


Fuil  comprehension  of  the  complete  program,  including  the  interaction  of 
the  several  subroutines,  can  be  obtained  only  from  a  detailed  study  of  the 
program  which  is  presented  in  its  entirety  in  the  appendix,  together  with  the 
discussions  of  the  several  subroutines  that  are  contained  in  the  sections 
that  follow  herein.  However,  it  is  informative  to  note  here  the  purposes  of 
each  of  the  several  subroutines.  These  purposes,  briefly  stated  and  without 
discussion  and  explanation  are  as  follows: 

INIT  "  Initializes  the  variables  to  their  proper  initial 

values. 

LOCATE  -  Locates  the  blast  pressure  pulse  in  relation  to  the 

masses  of  the  model. 


Performs  the  numerical  integration  of  the  equations 
of  motion  to  yield  changes  in  the  mass  point  coor¬ 
dinates  during  the  time  interval,  At,  resulting 
from  the  application  of  the  blast  forces. 


DEFORM 

EXTERN 

EXCITE 

YIELD 


Computes  the  internal  resistances  of  the  structure 
that  are  consistent  with  the  mass  point  displace¬ 
ments  that  are  computed  in  BETA. 

An  adjunct  to  DEFORM  to  assure  compatab 1 1 i ty  of  the 
axial  forces  and  moments  computed  in  DEFORM. 

Computes  the  soil  spring  forces  that  are  generated 
by  the  mass  point  displacements  determined  in  BETA. 

Maintains  the  time  histories  of  the  strains  in  the 
concrete  and  steel  in  each  of  the  several  elements 
of  the  model. 


LHtu‘  ‘  Determines,  on  the  basis  of  criteria  input  to  it, 

whether  failure  has  occured. 

As  noted  earlier,  a  complete  printout  of  the  entire  program,  including 
the  eight  subroutines,  is  presented  in  the  appendix.  A  list  of  the  FORTRAN 
variables  used  in  the  program,  along  with  their  definitions,  Is  also  given. 

For  simplicity  of  presentation  and  discussion,  the  variables  used  in  the 
following  discussions  of  the  several  subroutines  are  frequently  different  from 
their  FORTRAN  counterparts.  Therefore,  the  variables  are  defined  when  first  used 
in  program  development  and  are  summarized  for  convenience  in  alphabetical  order 
in  the  front  of  this  report. 
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2.  Subroutine  BETA: 

The  subrottine  BETA  performs  the  numerical  integration  required  by  the 
equations  of  motion  in  the  dynamic  portion  of  the  analysis.  The  method  used 
is  the  "Beta  Method"  developed  by  Newmark^  at  the  University  of  Illinois. 

Initial  conditions  of  displacement,  velocity,  and  acceleration  are  assumed 
to  be  known  at  time  t.  The  acceleration  of  the  mass  fit  time  (t+At)  is  then 
assumed,  and  the  corresponding  velocity  and  displacement  of  the  mass  at  time 
(t+At)  are  computed.  The  soil  and  structural  resisting  forces  that  correspond 
to  this  computed  displacement  are  then  determined,  which  when  used  together 
with  the  applied  force  at  t+At  make  it  possible  to  compute  an  acceleration  at 
that  time.  If  this  computed  acceleration  at  t+At  is  in  agreement  with  the 
assumed  acceleration,  the  assumed  acceleration  is  correct  and  the  resulting 
structural  configuration  is  the  actual  configuration  at  time  t+At;  If  not, 
the  computation  cycle  must  be  repeated  until  agreement  is  obtained.  The  dis¬ 
placement,  velocity,  and  acceleration  thus  determined  at  t+At  are  then  taken  as 
a  new  set  of  initial  conditions  and  the  procedure  is  used  to  determine  the 
conditions  at  the  end  of  the  next  time  interval. 


The  equations  used  in  the  Beta  Method  are: 


.(t+At)  .  *t  +i  (M)J;t  +J  (At)X(t+At)  (1) 

9(t+At)  .  ^  +  J  (At)v-t  +i  (At)y(t+At)  (2) 

Ax ft+At)  -  (At)  x  j  +  ^  -  0)  (At)2  xf  +  0(At)2  x{t+A°  (3) 

Ay[t+At)  -  (At)*{  +  -  0)  (At)2  y|  +  0(At)2  y|t+M)  (4) 

in  which 

AX(t*-At)  an(J  (t_f4t)  are  displacements  of  mass  i  in  the  x  and  y 

1  '  directions  during  the  time  interval  At. 

i" L"" ' r 

Newmark,  N.  M.,  "A  Method  of  Computation  for  Structural  Dynamics,"  Trans¬ 
actions  American  Society  of  Civil  Engineers,  Vol.  127,  1362 . 
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are  velocities  of  mass  i  in  the  x  and  y 
directions  at  the  beginning  cf  the  time 
interval,  and  are  known  either  as  initial 
conditions  or  from  analysis  of  the  preceding 
time  interval. 


x 


(t) 


and 


are  accelerations  of  mass  I  *n  the  x  and  y 
directions  at  the  beginning  of  the  time 
interval,  and  are  known  either  as  initial 
conditions,  or  from  computations  of  the 
preceding  time  interval. 


3 


is  the  dimensionless  constant  used  in 
numerical  Integration. 


The  value  of  the  dimensionless  parameter  3  depends  on  the  manner  in  which  the 
acceleration  is  assumed  to  vary  between  time  t  and  t+At.  In  the  problems 
solved  in  connection  with  this  report  3  *  1/6  was  used.  This  corresponds  ^o 
a  linear  variation  of  acceleration  within  the  time  interval.  The  time  Incre^ 
ment  At  must  be  chosen  carefully  in  order  to  insure  stability  and  convergence 
of  the  numerical  process. 


In  the  process  of  solving  the  equations  of  motions  listed  above  it  is 
clearly  necessary  that  the  blast  forces,  soil  spring  forces,  and  internal 
structural  resistances  be  known  in  order  to  compute  the  several  accelerations 
that  are  needed.  These  forces  are  computed  In  subroutines  DEFORM  and  EXCITE 
which  are  discussed  in  detail  later. 

In  programming  the  Beta  Method  for  this  problem  a  variable  time  interval 
At  was  used.  The  number  of  cycles  required  for  convergence  is  used  as  the 
criterion  for  the  choice  of  At.  If  the  problem  does  not  converge  in  eight 
cycles,  the  time  interval  is  reduced  to  half  its  value  and  cycle  is  restarted 
at  time  t»t;  on  the  other  hand,  if  convergence  is  reached  in  less  than  four 
cycles,  the  time  interval  is  doubled  for  the  next  time  increment.  In  this 
manner,  both  convergence  and  stability  of  the  method  are  assured. 


Convergence  in  both  the  x  and  y  directions  at  all  mass  points  in  the 
model  is  necessary  for  the  determination  of  the  structural  configuration  at 
time  t+At.  Absolute  convergence  is,  of  course,  not  required;  it  is  necessary 
only  that  the  acceleration  computed  to  exist  at  time  t+At  be  sufficiently 
close  to  that  which  was  assumed  to  maintain  acceptable  precision  in  the 
analysis. 
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Therefore,  convergence  can  be  checked  by  considering  the  difference  between 
the  assumed  and  calculated  accelerations  at  t+At,  without  regard  to  sign,  as 
fol lows: 


^  *i  ^  *1  Assumed  *i  Calculated  ^ 


where  1.  is  either  x,  or  y.. 
i  1  i 

If  Az,  is  less  than  a  prescribed  percentage  of  z  Qa icuia ted  t*ien 

acceleration  of  mass  i  In  the  direction  of  interest  is  considered  to  have 
converged.  This  check  is  applied  tc  all  masses  in  both  the  x  and  the  y 
directions.  If,  however,  Az.  is  greater  than  the  permissible  difference 

expressed  as  a  percentage  of  z . ,  it  is  still  possible,  if  z.  is  very  small, 

that  satisfactory  convegence  may  have  been  reached.  This  condition  is 
easily  checked  by  specifying  an  upper  acceptable  limit  for  the  absolute  value 
of  Az..  If  A*z'j  for  all  masses  is  smaller  than  either  of  these  allowable 

maximum  differences,  convergence  is  assumed  to  have  occurred,  and  the  analysis 
progresses  to  the  next  time  interval.  On  the  other  hand,  if  any  of  the 
accelerations  in  the  x  or  y  directions  of  any  of  the  masses  fail  both  the 
relative  and  absolute  checks,  new  accelerations  are  assumed  in  the  x  and  y 
directions  at  each  of  the  masses  and  the  analysis  is  repeated. 

Clearly,  the  values  of  Az.  that  are  considered  acceptable  depend  upon 
the  accuracy  sought.  Although  easily  changed,  the  program  as  currently 
written,  specifies  that  Az’j  shall  be  equal  to  or  less  than  0.1  in/sec/sec, 

or  0.05  percent  of  z’j,  whichever  is  larger. 


3.  Subroutine  CHECK: 

CHECK  is  a  subroutine  which  is  used  to  determine  whether  the  structure 
has  failed.  It  does  this  by  checking  the  strain  intensities  in  the  concrete 
fibers  that  are  used  to  approximate  the  cross  section  of  the  structural  element, 
and  it  is  called  after  acceleration  convergence  has  been  obtained  on  all  masses 
in  each  time  interval. 

Conceptually,  structural  collapse  is  assumed  to  have  occurred  when  any 
one  of  the  following  conditions  exist. 

1)  A  sufficient  number  of  yield  hinges  form  in  the  structure  to  trans¬ 
form  it  Into  a  mechanism 

2)  Crushing  of  concrete  occurs  in  any  fiber  at  the  number  of  locations 
defined  in  (1)  above. 

3)  A  combination  of  yield  hinges  and  concrete  fiber  crushing  exists  at 
the  number  of  locations  defined  in  (l)  above. 


When  the  criteria  indicate  that  the  structure  has  failed,  the  computer  prints 
where  and  how  it  failed  and  calls  an  error  routine  from  the  systems  library  to 
shut  off  the  computer.  As  presently  written,  the  program  assumes  a  two-hinged 
structure;  hence,  the  number  of  yield  hinge  or  crushing  locations  deflninq 
failure  is  two.  a 


4.  Subroutine  DEFORM: 

DEFORM  is  a  subroutine  which  computes  axial  bar  forces,  and  moments  and 
shears  at  the  joints  from  changes  in  the  mass  point  coordinates  and  the  stress- 
strain  curves  for  concrete  and  steel.  The  subroutine  uses  coordinate  changes 
to  compute  strains,  which  are  then  converted  to  stresses  and  finally  to  forces, 
moments  and  shears.  It  is  patterned  closely  after  a  similar  program  developed 
by  Dawkins.^ 

Figure  3  shows  bar  i  connecting  joints  i  and  (l+l)  at  two  Instants  of 
time.  During  the  time  interval  At  the  joints  have  moved  to  new  positions 

causing  coordinate  changes  of  Ax!+At,  Ax£f,  Ayf*At,  8nd  AyfJJ*.  The  coor¬ 
dinates  of  joints  i  and  (i+1)  at  time  t+At  are,  therefore: 

-  xf  +  Ax™;  y™  .  yt  +  Ayt «4t 


+  Ax™*  v™  -  +  Av™ 

*1+1  *  y|+l  y|+l  +  Ayi+1 


Using  these  new  coordinates  the  slope  of  bar  i  at  time  t+At,  a!™  is  defined 
bV  ..t+At  t+At  ’ 


-  ARC7AN 


M+l  7i 
t+At  t+At 
*1+1  "  xl 


Then,  the  incremental  change  in  ff(,  Ao™,  during  the  interval  At,  Is 


puted  by 


<r4t  -  [  «4t  -  4«:+4t 


SIN  » 


- «!'  -  )  cos  <l|  ♦ 4L;+“)  0 

The  incremental  change  in  the  length  of  bar  i,  AL™,  during  the  Interval  At, 
is  computed  from  1 


Dawkins,  W.  P.,  "Dynamic  Response  of  a  Tunnel  Liner-Packing  System",  Ph  D 
Thesis,  University  of  Illinois,  1966.  9  V  *  '  ' 
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At  t  =  t 


t  8  t  ♦  At 


*/ 


W  I 


FIGURE  3 .  DISPLACEMENT  OF  BAR  DURING  TIME 
INCREMENT 


AL™'  -  (>x£f  -  Ax™)  COScv™ 

*■  (Ayi+i 1  ■  AyrAt)  sin  »rAt  (9) 

It  Is  Important  to  note  that  the  equations  above  give  values  for  changes 
which  occur  during  the  time  Interval  between  t  and  t+At.  To  obtain  the  total 

TOTAL 

change  in  bar  length,  ALj  It  is  necessary  only  to  sum  the  previously  com¬ 
puted  incremental  values  of  ALj  as  follows: 
t"t+At 

TOTAL  _  <->  „n 
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» 


Similarly,  the  total  rotation  of  bar  i ,  Aa,  ,  is  obtained  as  the  sum  of 
the  incremental  rotations. 


t”t+At 

to!™1 .  I  to,  do 

t*0 

From  these  values  of  total  bar  rota  Jon  and  change  in  length,  strains  and 
curvatures  can  be  computed.  The  average  strain  in  bar  I  is 


AVG 

BARI 


»  AL¬ 


TO  TAL 

i 


(12) 


where  L  is  the  original  bar  length  and  carries  no  subscript  since  all  the 
bars  originally  have  the  same  length. 


It  is  assumed  initially  that  all  rotations  are  concentrated  at  the  Joints; 
that  is,  the  bars  rotate  as  rigid  bodies.  Figure  4  shows  Joint  I  with  portions 
of  bars  i-1  and  i  attached.  Bar  i  has  undergone  a  positive  angle  change 

(clockwise  rotation)  AarJ°TA1,  while  bar  1-1  shows  a  negative  angle  change 

(counter-clockwise)  AqJ®Ta*-.  The  rotations  of  these  two  bars  results  in  a 

*■  1 " 1  t 

decrease  p~,  the  angle  between  the  two  bars.  This  angle  change,  Ap|>  can  be 

computed  from 


=Aq 


TOTAL 

i 


-  ta 


TOTAL 

i-1 


(13) 


equation  (13)  gives  a  positive  value  to  Ap*  when  the  intersection  angle  is 
reduced.  An  average  value  of  curvature  at  the  Joint,  l\T,,  is  obtained  by 

4  t  ^  * 

assuming  the  total  angle  change,  Ap.,  to  be  uniformly  distributed  over  a 
distance  equal  to  one  bar  length.  1 


♦jT!  *  *pi/L  (14) 

Equations  (13)  and  (14)  are  consistent  with  a  sign  convention  for  which  a 
decrease  in  the  central  angle  pj  produces  a  positive  curvature  at  the  joints. 

By  defining  a  positive  moment  as  one  which  produces  compressive  strains  on 
the  outside  of  the  structure,  moments  in  the  model  are  consistent  with  those 
given  by  the  differential  equation  for  bending, 

M  «  _  ei  $-1 

dx 


19 


FIGURE  4.  RELATIVE  BAR  ROTATIONS  AT  MASS 


where  M  is  the  bending  moment,  E  the  modulus  of  elasticity  of  the  material 

comprising  the  section,  I  the  moment  of  inertia,  and  the  curvature  of  the 
section.  dx^ 

H  .  :  "  M‘ 

The  preceding  equations  yield  average  axial  strains  in  the  bars  and 
average  curvatures  at  the  joints  of  the  model.  However,  because  of  the  nature 
of  the  model  and  the  fact  that  inelastic  response  of  the  system  Is  to  be 
considered,  it  was  necessary  to  develop  a  procedure  to  account  for  the  inter¬ 
action  of  axial  force  and  moment  both  at  the  joints  and  in  the  bars.  Thus, 
the  axial  force  in  bar  i  was  determined  from  the  average  axial  strain  computed 
from  Equation  (12)  for  the  bar  together  with  bending  strains  Which  were  assumed 
to  be  the  average  of  the  bendinq  strains  defined  by  the  curvatures  at  Joints 
i  and  i+1  as  given  by  Equation  (14).  Similarly,  the  moment  at  joint  i  was 
computed  1  *om  the  strains  given  by  the  curvature  at  that  Joint,  Equations  (14), 
together  w  th  an  axial  strain  which  was  assumed  to  be  the  average  of  the  axial 
strains  ot  the  two  adjacent  bars  as  given  by  Equation  (12).  The  following 
equations  are  consistent  with  these  assumptions. 
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The  average  axial  strain  in  joint  i,  €.—  »  Is  taken  as 

1 1 

AVG  AVG 
AVG  *BAR1  CBAR  1-1 
*JTi  2 


and  the  curvature  at  mid- length  of  bar  i,  ♦bar i *  *s  assume^  to  be 

♦jTi  +  *JT  i+1 
b.'Ri  2 


Assuming  the  angle  chsnges  to  be  sufficiently  small  that  the  sine  of  the 
angle  Is  approximately  equal  to  the  angle  expressed  in  radians*  and  adopting 
a  sign  convention  that  is  positive  when  z,  a  distance,  is  above  the  plastic 
centroid  the  effects  of  axial  and  bendfng  strains  can  be  combined  using  the 
following  expressions: 


z  AVG 

*BARi  "  cBARi  -  2  ▼nARI 


JTi  JTi 


*.it  i  i  2  * 


Z  Z 

where  «BARi  and  CjTJ  are  strains  in  bar  I  and  joint  I  at  some  distance  z  above 

or  below  the  plastic  centroid.  Figure  5  shows  this  strain  summation  for  bar  I. 
For  numerical  couputat ion,  the  bars  and  joints,  whose  cross  section  is  shown 
in  Figure  6,  are  represented  by  a  number  of  fibers  as  shown  In  Figure  7.  The 
(•esulting  equations  for  the  strains  in  fiber  j  of  joint  i  and  bar  I  are 


AVG  _  Ijj  /IhJjJJ  ]  . 

C  BARi  L2  2  J  nfJ  ▼£ 


€AVG  .  [H  .  fUblx  H  1 
JTI  (.2  V  2j  nfJ  *JYi 


in  which  eBARj,  and  e^j.  are  the  strains  in  fiber  j  of  bar  i  and  joint  i,  H 

is  the  total  depth  of  the  cross  section  and  n,  is  the  number  of  fibers  used  to 
approximate  the  cross  section.  ' 
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FIGURE  5.  STRAIN  SUMMATION  IN  BAR 


The  stress  distribution  can  now  be  obtained  from  applicable  stress- 
strain  curves.  Figure  8  is  an  idealized  trllinear  stress-scraln  curve  of 
concrete,  on  which  the  following  notation  is  used. 

E  -  modulus  of  elasticity 

v 

*cr  *  strain  at  which  concrete  begins  to  fail 

«u  *■  ultimate  strain  or  strain  at  which  the  concrete  crushes 

*o  »  strain  axis  Intercept 

The  capitalized  names  in  parentheses  are  the  FORTRAN  names  assigned  to  the 

variables.  The  concrete  is  assumed  to  unload  along  a  line  parallel  to  £  , 

c 

its  initial  modulus  of  elasticity.  The  tensile  strength  of  concrete  is 
neglected. 
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FIGURE  7.  FIBER  STRAIN  REPRESENTATION  FOR 
COMPUTER  ANALYSIS 
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The  axial  forces  in  the  bars  and  joints  are  then  computed  from  the 
stresses  using  the  following  equations: 


■  1  (oJ8AR  •  V 


F0RCEJT|C  ■  l  Mt,  •  V 

J-l 


where 


FORCeE?b^»  F0RCE*i~C  are  the  forces  developed  in  the  concrete  of  bar 
BARI  JT1  and  joint  i 


J  J 

^BARi*  °JTI 


are  the  stresses  in  fiber  J  of  bar  I  and  joint  t, 


Similarly  the  moments  in  the  Joints  and  bars  are  computed  from 


^c-i{R,-p-e¥)y}  o» 

J-1 

in  which  and  MOM*'.??*'  are  the  bending  moments  developed  by  the  concrete 

BARI  Jll 

bar  in  bar  i  and  joint  I . 

The  preceding  analysis  has  neglected  the  reinforcing  steel  in  the  bars 
and  joints.  The  steel  is  distributed  equally  in  two  layers,  one  near  the 
inner  and  the  other  near  the  outer  face  of  the  cross  section;  It  i$  placed 
symmetrical ly  about  the  plastic  centroid.  It  is  assumed  to  behave  as  a  purely 
elasto-plastic  material  shown  In  Figure  9. 
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On  the  basis  of  the  preceding  analysis  of  concrete  strains,  the  strains 
in  the  reinforcing  steel  of  bar  i  are 


OUTER  STEEL 


«AVG  +il.  * 

*BAR i  2  T8ARi 


INNER  STEEL  o 


AVG 

*BARi  “  2  fBARi 


in  which  d1  is  the  distance  between  the  layers  of  reinforcement.  Similar 
equations  exist  for  the  steel  in  Joint  i. 

Stresses  in  the  steel  are  determined  from  the  stress-strain  diagram  of 
Figure  9.  Forces  in  the  reinforcement  are  computed  by  multiplying  the  area 
of  the  reinforcement  by  the  stresses.  If  any  of  the  steel  is  in  compression, 
it  Is  necessary  to  use  an  equivalent  area  which  is  given  by 

A  „  (n-.O.  A  (28] 

AEQUIV  n  s 

in  which  n  is  the  ratio  of  Ec  to  E  ,  A  is  the  area  of  the  reinforcing  steel 

S  C  5 

md  A-n||TU  is  the  modified  compression  area  of  the  reinforcement. 
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Then,  combining  the  effects  of  concrete  and  steel,  the  total  force  in 
bar  I  is  computed  by 

FORCE  *»  F0RCEE^G  +  FORCE^^"^  +  F0RCE^^ER  ( 

F0RCEBAR!  F0RCEBARi  F0RCEBARI  F0RCEBAR1  ' 


and  the  total  moment  in  joint  i  Is 


d*  Q 


FORCE 


INNER  STL 


-  FORCE 


OUTER  STL" 


The  transverse  shear  in  bar  i,  which  corresponds  to  the  difference 
between  the  moments  at  its  ends,  is  then  computed  as 


’JT  i+1 
L, 


, -pupii'  i  ^  n  iiuimr  r  fl  1 


Figure  10  shows  bar  i  with  positive  thrusts,  shears  and  moments  acting 
upon  it.  These  are  the  values  obtained  from  equations  (29),  (30),  and  (31), 

It  is  of  interest  to  note  that  the  method  presented  above  cannot  always 
be  relied  upon  to  locate  the  neutral  axis  of  the  cross  section  properly.  As 
an  example,  consider  a  reinforced  concrete  beam  subjected  to  pure  bending. 
Using  the  aboveanalys is,  a  moment  and  a  compressive  thrust  would  be  computed. 
T  is  thrust,  which  in  actuality  is  non-existent,  arises  because  the  analysis 
places  the  neutral  axis  at  the  plastic  centroid.  The  concrete  in  tension 
produces  no  force  while  the  compressed  concrete  develops  a  force.  The  re¬ 
inforcement  develops  both  tensile  and  compressive  forces,  but  an  unbalanced 
compressive  force  remains  on  the  net  cross  section.  By  adjusting  the  strains 
until  they  produce  no  thrust  the  neutral  axis  could  be  properly  located. 

EXTERN,  a  subroutine  described  later,  calculates  a  thrust  in  the  bars  and 

joints  directly  from  the  average  strains  and  These  arc  correct 

values  and  the  strain  distributions  in  the  method  described  above  are  adjusted 
until  the  thrusts  computed  in  this  manner  agree  within  acceptable  limits,  with 
those  given  by  EXTERN.  The  corresponding  moments  are  then  assumed  also  to  be 
correct.  Returning  to  the  example  of  the  beam  subjected  to  pure  bending, 

EXTERN  would  calculate  a  thrust  equal  to  zero  since  the  ends  of  the  beam  have 
not  been  displaced  in  an  axial  direction.  DEFORM  would  then  use  the  values 
calculated  by  EXTERN  to  modify  the  fiber  strains  until  the  computed  thrusts 
are  within  acceptable  limits,  equal  to  zero.  This  correction  is  required  only 
if  the  strain  changes  sign  between  the  outer  and  inner  faces  of  the  cross 

section.  The  problem  does  rot  exist  if  the  entire  section  is  in  either  tension 
or  compression. 


5.  Subroutine  EXCITE: 

The  blast  and  soil  spring  forces  acting  upon  the  structure  are  computed 
by  EXCITE.  Deformations  in  the  springs  are  used  to  determine  spring  fo,  ces 
and  a  pressure-time  diagram  is  used  to  determine  blast  pressures. 

LOCATE,  a  subroutine  discussed  later,  calculates  tht  total  time  that 
has  elapsed  since  the  wave  first  engulfed  each  of  the  masses  on  the  structure. 
Using  a  pressure-time  diagram  such  as  Figure  11  and  the  elapsed  time  since 
initial  contact  of  the  mass  by  the  blast  wave,  it  is  possible  to  determine  the 
pressure  level  acting  on  mass  i.  The  pressure  from  the  diagram  is  the  vertical 
pressure  acting  downward  on  mass  i  at  the  time  in  question.  The  pressure  on 
mass  i  is  assumed  to  extend  over  half  a  bar  length  on  each  side  of  it.  To 
obtain  the  vertical  blast  force  on  mass  i  the  pressure  on  that  mass  is  multiplied 
by  the  horizontal  projections  of  half  of  bars  i-1  and  i.  Therefore, 
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”her«  FO«CE®^ST  is  the  concentrated  bias,  force  action  a 

and  PRESSURE,  Is  the  vertical  blast  Awards  on  mass  1 

blast  force  „  ...  :  .  "  'C""*  ““  *  *•  <«• 
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horizontal  free-field  pressures  ^hTch ' al  ly  the  same  manner  from  the 
constant  (RATIO)  times  the  cor»2-  K  h  ^  ?SSUmed  to  be  equal  to  a 
the  horizontal  blast  force  on  mass  I  !"9  }'ert,cal  free-field  pressures.  Thus, 
projections  of  half  of  bars  i-1  and  f  *  preSSure  mult,Pn«d  by  the  vertical 
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soil.  0  of  hor,z°"tal  to  vertical  blast  pressures  in  the 


using  the^fol lowing  'varfables'!''  65  *"  Fl9UrC  Is  read  ,nto  th«  computer 


A1  "  first  pressure  amplitude 

A2  -  second  pressure  amplitude 

TR1  -  time  at  which  pressure  attains  amplitude  A1 

TAR2  -  time  at  which  second  pressure  wave,  A2,  arrives 

TR2  -  time  at  which  pressure  attains  amplitude  A2. 


If  A1  A2  a  step  pulse  of  infinite  duration  result*  Tf  tabs  tb-» 

line  connecting  A1  and  A2  results.  These  var ion! 'J 7  ™2,  *  $traI9ht 

pulse  are  shown  in  Figure  12.  Necessarv  f  thf  bas  ,c  Pr«S5'jr« 

TR2  are  that  Necessary  restrictions  on  choice  of  TRI,  TAR2  and 


TR1  <  TAR2 
1AR2  <  TR2. 


The  resistance  of  the  soil  surrounding  the  structure  is  aDoroximated  h« 
springs  attached  to  the  messes.  Since  mos?  soP  hes  UUelfLv  ^.lf. 

llrulr a”^d  t0  «"**  I"  compression!^  As^noted "* 
anticipation  “  ! TtL'l  b'.: ""“fir  hoover, 


~  resistance  is  assumed  to  be  linearlv  elastlr 

anticipation  ->f  extending  the  program  to  treat  multilinear  soil"  resistance 

^"fillS.  Bltial  eUStiC  SlOPe  iS  defined  b*  tha  sPr^ng°st Iffness*nCe 
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FIGURE  12.  VARIATIONS  OF  BLAST  LOAD  FUNCTIONS 


where  FORCYY(l)  Is  the  yield  pressure  of  the  springs  at  mass  1  and  DEFYY(I) 

Is  the  yield  deformation  of  the  springs  at  mass  i.  Since  an  elastic  spring 
is  being  used,  DEFYY(I)  has  no  physical  meaning  and  Is  used  only  to  define 
k  .  If  DEFYY(I)  -  1.  is  read  into  the  computer,  then  k.  -  FORCYY(l).  Figure 
13  is  the  pressure-deformation  diagram  of  the  springs  used  to  approximate  the 
soi  1 . 

The  vertical  force  on  mass  i  from  the  vertical  soil  spring  I? 
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and  the  horizontal  force  on  the  mass  from  the  horizontal  spring  is 


FORCE 


SPRING 


-  k  Ax.  (---  SIN  a._,  +  ^  SIN  or,) 


The  values  of  Ay.  and  Ax.  are  the  total  deformations  in  the  springs  and  are 
equal  to  the  total  coordinate  change  of  mass  i  as  calculated  by  BETA. 


Prttsur* 

<  poi ) 


FORCYY(I) 


DEFYY(I) 


Deformation  (la) 


FIGURE  13.  SOIL  RESISTANCE  FUNCTION 


6 .  Subroutine  EXTERN : 


The  purpose  of  EXTERN  has  already  been  discussed  In  some  detail  In  the 
section  on  DEFORM.  In  EXTERN,  the  average  bar  and  joint  strains  are  used  to 
calculate  thrusts  which  are  then  used  In  DEFORM  to  correctly  locate  the  neutral 
axis  of  the  section,  and  thence,  the  bending  moment  on  the  section. 


The  axial  force  in  the  concrete  of  bar  I  Is  given  by 


F0RCEC0MC  -  oCa0NC  AVG 
BARI  BARI 


where  H  is  the  section  depth  and  crG!?^  is  th*  concrete  stress  which  Is 
obtained  from  Figure  8  to  be  consistent  with  «R«nj  from  Equation  (12). 


Similarly,  axial  force  in  the  steel  reinforcement  Is 
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WAVG*  CT|3^Ri  iS  thC  Stee*  stress  wh^ch  Is  obtained  from  Figure  9  using 

*BARI  from  Equation  (12).  If  the  bar  is  In  compress  ion  A  ,  must  be 
modified  by  Equation  (28). 

The  total  axial  force  In  bar  i  is  the  sum  of  the  axial  forces  carried  bv 
the  concrete  and  steel,  or 
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(39) 


Similar  equations  can  also  be  developed  for  the  axial  forces 


in  Joint  i . 


7 .  Subroutine  INIT: 

The  variables  in  the  program  are  initialized  to  their  proper  values  by 
INIT.  The  subroutine  is  called  only  once  at  the  beginning  of  the  program. 


8.  Subroutine  LOCATE: 


The  purpose  of  this  subroutine  is  to  locate  the  advancing  pressure  pulse 
with  respect  to  the  model  under  study.  Two  versions  of  LOCATE  were  developed. 
One  version  was  developed  to  study  the  special  case  of  a  horizontal  wave, 
propagating  vertically  downward  as  shown  in  Figure  14;  the  other  version  con¬ 
sidered  the  more  general  case  of  an  inclined  wave  propagating  horizontally 
across  the  structure  as  shown  in  Figure  15.  In  its  limit,  the  latter  case 
can  be  made  to  approximate  the  former,  which  is  a  special  case  developed 
primarily  to  facilitate  the  checking  of  the  computer  program.  A  symmetrical 
structural  system  subjected  to  the  wave  generated  by  the  vertical  version  of 
LOCATE  responds  symmetrically;  this  fact  substantially  simplified  the  pro¬ 
gram  checking  studies. 

If  a  study  of  structural  response  to  unsymmetr ical  loading  is  of  in¬ 
terest,  the  second  version  of  LOCATE  should  be  used.  In  any  case,  only  one 

of  the  versions  of  this  subroutine  should  be  compiled  with  the  other  sub¬ 
routines  to  build  a  complete  program. 

Both  versions  of  LOCATE  use  the  coordinate  of  mass  i  to  determine  the 
length  of  time  that  the  pressure  has  been  acting  upon  that  mass.  The  vertical 
version  of  LOCATE  uses  y.  of  mass  i  for  location  purposes,  while  the 
horizontal  version  uses  x..  The  location  method  for  the  horizontal  version 
is  similar  to  that  of  the  vertical  and  will  not  be  discussed  here. 

During  a  time  interval  At,  a  wave  propagating  downward  across  the  struc¬ 

ture  moves  a  distance  Ay  which  Is  determined  by 

Ay  =  VERTICAL  VELOCITY  x  At  (40) 

where  the  vertical  velocity  of  the  wave  front  is  read  into  the  computer  as  part 
of  the  input  data.  The  wave  makes  contact  with  the  top  of  the  structure  at 
time  t«0;  thus,  at  this  instant,  the  wave  front  is  at  the  highest  mass  on  the 
structure,  and  is  designated  as 
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FIGURE  14.  BURIED  ARCH  UNDER  VERTICAL  SHOCK 
FRONT 


FIGURE  15.  BURIED  ARCH  UNDER  INCLINED  SHOCK 
FRONT 
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At  time  t*  QFAt,  the  location  of  the  front  is 


WAVE  FRONT WAVE  FRONT1*0  -  Ay 


or,  in  general,  at  any  time  t«t, 

rt«t 


WAVE  FRONT 


WAVE  FRONT1*1"*1-  Ay 


(42) 


To  determine  the  amount  of  time  that  the  wave  has  been  acting  on  mass  i, 
the  distance  between  mass  i  and  the  wave  front,  6yf*t,  must  be  determined. 
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(43) 


If  this  value  is  negative,  the  wave  front  is  above  masr,  i  and  the  computer 
value  of  T(l) ,  the  elapsed  time  since  wave  contact  with  the  mass,  is  set 

equal  to  zero.  If  fiyj”*  is  positive,  the  value  of  T(l)  is  obtained  by 


T(I)  -  6y;=t/  (VERTICAL  VELOCITY) 


(44) 


The  values  of  T(l)  computed  in  LOCATE  in  this  manner  are  used  in  EXCITE. 

The  pressures  on  the  several  masses  are  determined  at  times  T(l)  from  a 
pressure-t ime  diagram  that  is  stored  by  the  computer  in  that  subroutine. 

9.  Subroutine  YIELD: 

The  time  histories  of  the  strains  in  the  various  parts  of  the  structure 
are  computed  and  stored  in  YIELD.  The  histories  of  the  strains  in  both  the 
concrete  fibers  and  in  the  steel  reinforcement  are  necessary  because  neither 
of  the  materials  is  purely  elastic,  the  determination  of  stress  from  strain 
in  a  vibration  problem,  where  there  is  periodic  loading  and  unloading,  is 
impossible  without  a  complete  strain  history.  Such  histories  are  necessary 
so  that  the  accumulated  plastic  sets  can  be  properly  considered  when  stresses 
corresponding  to  computed  strains  are  read  from  stress-strain  diagrams  such 
as  those  shown  in  Figures  8  and  9. 

YIELD  checks  every  strain  in  the  concrete  fibers  and  the  reinforcing 
steel  of  the  several  bars  and  joints  of  the  model  to  determine  if  and  when 
the  strain  exceeds  the  yield  strain  for  the  material  concerned.  If  the  yield 
strain  has  been  exceeded  it  calculates  the  plastic  set,  or  permanent  deformation, 
by  subtracting  the  yield  strain  from  the  total  strain.  In  bar  I,  the  plastic 
set  in  fiber  j  of  the  concrete  is  given  by 
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Similarly,  in  the  outer  and  inner  reinforcing  steel  the  plastic  sets  are, 
respectively 
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CONC  STEEL 

where  ey  and  «  are  the  yield  strain*  In  concrete  and  steel. 


subroutine  checks  the  strains  before  BETA  starts  a  new  time  Incre¬ 
ment  and  modifies  the  plastic  sets  In  accordance  with  the  lost  computed  Itralhs 


SECTION  IV 


ILLUSTRATIVE  APPLICATION  OF  PROGRAM 


1 .  Description  of  Soil-Structure  System: 

To  illustrate  the  applicability  of  the  computer  program  developed  in 
this  study,  it  was  used  to  investigate  the  Influence  of  localized  variations 
in  soil  resistance  on  the  blast  pressure  that  is  required  to  damage  a  typical 
buried,  pin-ended,  reinforced  concrete  arch. 

The  arch  that  was  analyzed  is  illustrated  in  Figure  16(a),  It  had  a 
radius  of  15  ft.  and  a  total  thickness  of  13.5  inches,  which  corresponds  to  a 
thickness-radius  ratio  of  0.075  The  central  angle  of  the  arch  was  180°, 
and  it  was  reinforced  with  1.0  percent  steel  on  each  face.  The  steel  was 
placed  circumferentially,  1  1/2  inches  from  each  face  of  the  section.  A 
cross  section  through  the  arch  is  shown  in  Figure  16(b),  and  the  stress-strain 
diagrams  for  the  concrete  and  steel  were  idealized  as  shown  in  Figures  16(c) 
and  (d)  . 

For  purposes  of  analysis,  a  section  of  the  arch  1.0  inch  long  was  con¬ 
sidered.  This  soil-arch  system  was  represented,  as  shown  in  Figure  17(a), 
by  a  six-bar  system,  with  the  masses,  resistances,  and  forces  concentrated 
at  the  node  points.  Each  bar  was  93.175  inches  long,  and  had  the  same  cross 
section  as  did  the  arch  being  represented. 

All  soil  springs  were  assumed  to  be  linearly  elastic  and,  at  any  given 
mass  point,  both  the  horizontal  and  vertical  springs  were  assumed  to  have  the 
same  stiffness.  The  basic  soil  springs  stiffness  of  all  springs,  against  which 
the  comparative  effect  studies  were  made,  was  computed  as  described  in 
Section  III  to  be  consistent  with  a  foundation  modulus  of  50  ps!  per  inch 
of  deformation. 

The  masses  were  determined  on  the  assumption  that  the  unit  weight  of  the 
surrounding  soil  was  100  pcf,  and  that  soil  within  one  irch  radius  of  the 
arch  vibrated  with  the  arch.  Thus,  each  mass  was  computed  from  the  weight 
of  the  "effective11  adjacent  soil,  to  be  1.301  lb-sec2/in. 

The  system  was  subjected  to  a  pressure  pulse  of  the  form  shown  in 
Figure  17(b).  This  pressure  p-jlse  was  assumed  to  have  a  horizontal  front  and 
to  propagate  downward  through  che  soil  at  a  seismic  velocity  of  4000  fps. 
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Ill 

Ilf  ,(* 


(o)  Lumped  Mats  System 


Pressure  (psl) 


Alt  A2 


(b)  Blast  Pressure  Pulse 


Time  (sec.) 


FIG.  17.  SOIL” STRUCTURE  SYSTEM  IDEALIZED  FOR 
ANALYSIS  . 
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2 .  Influence  of  Variations  In  Sofl  Resistance: 

The  system  as  previously  described  was  first  analyzed  to  determine  the 
blast  pressure  level  that  was  required  to  produce  failure  under  a  uniform  soil 
resistance  corresponding  to  a  foundation  modulus  of  50  psl/in.  of  deformation, 
Failure  was  defined,  as  described  in  Section  III,  to  have  occurred  when  one  of 
the  following  conditions  existed. 

'I  -  .n',1 

(a)  Yield  hinges  developed  on  two  sections  in  the  model, 

(b)  Concrete  crushing  occurred  in  fibers  at  two  sections  in  the  model, 

(c)  A  yield  hinge  and  concrete  crushing  developed  at  the  same  or  at  two 
different  sections  iri  the  model. 


This  basic  blast  failure  pressure  was  found  to  be  95.5  psi. 

Analyses  were  then  carried  out  to  deterr.ine  the  extent  to  which  this 
blast  failure  pressure  would  be  reduced  because  of  equal  reductions  in  the 
soil  spring  stiffnesses  at  mass  points  2  and  f»,  the  other  soil  spring  stiff¬ 
nesses  remaining  unchanged.  The  results  of  these  analyses  are  shown  in 
Figure  l8(u)  . 

Similar  analyses  were  carried  out  to  determine  the  influence  on  the 
blast  failure  load  of  varying  the  soil  spring  stiffnesses  mass  points  3 
and  5,  the  other  stiffness  remaining  unchanged.  The  results  of  these  studies 
are  shown  in  Figure  18(b), 

The  results  portrayed  in  Figure  18  are  self-explanatory.  It  is  appro¬ 
priate,  however,  to  point  out  thut  these  data  are  indicative  of  the  kind  of 
information  that  can  be  obtained  through  the  use  of  the  computer  program  that 
was  developed  in  this  project.  If  properly  employed,  It  could  be  used  to 
yield  information  that  would  be  of  significant  value  in  the  evaluation  of 
hardness  of  existing  facilities  as  well  as  In  the  design  and  construction  of 
similar  facilities  In  the  future. 


Pm  *  Bloif  Prttturt  To  Product  Follurt 

P^s  Blast  Prttturt  To  Product  Folluro  Undtr 
Uniform  Soli  Rttitlonct 


Pm 


(a)  Influence  of  Variations  in  Soil  Resistances  at  Mass  Points 
2  and  6. 
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(b)  Influence  of  Variations  in  Soil  Resistances  at  Mass  Points 
3  and  5 . 

FIG.  18.  REDUCTION  IN  BLAST  FAILURE  PRESSURE  AS  A  FUNCTION 
OF  LOCALIZED  REDUCTIONS  IN  SOIL  RESISTANCE. 
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APPENDIX 


Th*  following  pages  contain  a  raproductlon  of  tha  actual  computar 
printout  of  th«  complete  program  that  was  developed  In  ti  Is  projact.  For 
ease  In  studying  It,  an  alphabetic  listing  of  all  FOR THAN  variables  used 
In  the  program,  together  with  the’r  definitions,  Is  presented  In  Sub¬ 
section  2.  Similarly,  the  basic  Input  data  are  summarized  In  Sub-section  3, 


1.  Computer  Program  Printout; 
a.  Ha  In  Program 

DIMENSION  PALFAI  10)*PX(11>.PY<11).X<11)«Y<’1)  ,DOARL<  101  *EPS0AW(  10) 

1  . TEPBAR  < 10)  i OALFA (10) , DAN  JT  { 1 1 ) .CRVJT( 1 1 ) *CRV0AR( 1C)* 

2  TF£PS8( 10*20)  .EPSJT( 1 1) *  TFEPS J (  1  1  *  20 )  . TS 

3EPBT (10),  TSEPBS(iO)  *TSEPJT(11>  ,TS 

4EPJBI 1 1 )  . PL SETS ( 10.20) «FS 1GB ( 10.20) 

OD I  MENS l ON  PLSETJ ( 1 l , 20 )  .  F  S I QJ (11*20)  .CFORCBt 10) *CMOM®< 10) .CPOMCJl 1 
1 1 )  « CMOM J ( 1 1 )  «  SEPYPT ( 1 1 >  .SEPYPCI 1  i > .ESEPBTUO) .SPSTgTi 10) *SFO«BT ( 10 
2) *SPSTBB( 10) ♦ ESEPBB ( 10) . SFORBB < 1 0 ) » SPST JT ( l l > ,ESEPJT( 1 1 ) **rORJT( 1 l 

3  )  ,  ESEP JB  <  1  1  )  *SPSTJ6<  1  1  )  ,SFO«JB< 1  1  )  .FORCBC  1  0  I  *  I3MOMB  (  1  0  )  *FO«C  J  (  11  >  *B 

4  MOMJ  (  11)  *V0(  10)  «  WFRTX  (  1  1  )  .  XDIFFXU  1  ) 

OD I  MENS  I  ON  FORCY (11)  *DEFYY(  1  1  )  •  DEFY ( 1 1 )  . DEFX 1 ( 1 1 ) *OEFX2<  1 l 

1  ) 

2* A (  1 1 )  *FY< 1 l  )  *FX1 < 1 1 ) *FX2<  1 1 )  *T< 1 1  >  «  EPSETB  t 10 ) *EFORCB< 10 ) .TEPJT( 1  I 
3  > .EPSETJ ( 1 1  )  .EFORCJ (11)  .ESPSTB (10) *EESEPB( 10) *EF0RS8< 10) .ESPSTJ) 1 l 
4)  • EESEP J ( 1 1 )  *EFORSJ( 1 l )  .COUNTB (10) . FORDFB  < 10) .D1FEPBC 10) 

DIMENSION  CMASS  < 10) *FV( 10) *FH( 10) *  VV  <  10) *  VH  ( 10) »FV1 ( 10) *FV 

12( 1 0) *FH1 ( 10) *FH2< 1 0) . VV1 ( 10) .  VV2( 10) *VMl < 10) . Vh2 (20) *DX( 10) »DY( 10 
2) *  DDX  <  10) *ODY( 10) *  ALFA ( 10) *0DXEND( 10)  ♦ DDYEND ( 10) 

3  *  DXEND  < 10) *  DYE NO (  10)  *F0RDFJ(  1 1 ) .DIFEPJl  1  1 )  (BLAST  1  (11) .BLAST 

42(11)  .  BLAST Y( 1 1 ) »  PDOX  <  10)  ♦ PDDY ( 1 0 ) »LC ( l 0 >  « JC ( 1 1 1  * LY H 0 ) « J Y < l 1 > 
DIMENSION  CHKYl (10) *CHKY2(  10) .XORIGI  l  1  )  *Y0RIG( 1 1  ) 

DIMENSION  DELX(IO) « DELY (10). TDELX ( 10>*T0ELY( 10) 

C 

OCOMMON  PALFA  *PX*  PY . X, Y  *  DBARL .EPSBAR ♦ TEPBAR. DALFA . OANJT.CRVJT .CRV8A 

IR.  TFEPSB.EPSJT,  TFEPSJ.  TSEP8T  , TSEPBB 

2  . TSEPJT  .TSEPJB  . PLSETB .FS I Gb  . PLSETJ . FS | GJ ♦ 

3C-ORCB.CMOMB.CFORCJ, CMOM J, SEPYPT. SEPYPC.ESEPoT.SPsTBT.SFORBT.SPSTB 
4B.  ESEPBB. SFORBB. SPSTJT  .ESEPJT ,  SFQR JT  .ESEPJB . »PST  JB  .  SFORJB *PORCe> ,BM 
50MB  «FOWC J  *BMOM J  »N  *  NN. BARL  *  NF 1 8RE . D . DPR  I M *  AS. EC . CEPSY . CEPSCR . CEPSO . 
6CEPS0.ES. SEPSy  . VB.WFRTX.XDIFFX 

OCOMMON  A1  .A2.TR1  .’•AR2.TR2.PHEE  .FORCY  .  DEF  YY  .  DEFY  .  Dc.Fx  1  .D 

1EFX2  .A.FY.FX1  .PX2*T.oUV£ 

2L. Q.DELT1M.EPSETB.EF0RCB.TEPJT. EPSETJ. EFORCJ. ESPSTB. EESEPB.EFORSU. 
3ESPST J, EESEP J .LFORSJ. COUNTB. FORDFB. D I FEPB .FOHOF J . D I FEPJ 
COMMON  HGT. CMASS. FV.FH.VV.VH.FV1 . FV2 , FH 1 .FH2.VV1 .VV2.VH1 .VH^. GAMMA 

I S .  DX . DY , DUX , UUY . AU  F A , DUXENO . DDYEND  • uxEND . DYLND . b • AH-nB . 

2  T I  ME. HAT  10. ALERR.ULAST1 .DLAST2 .BLAST Y 

C  OMMON  PDDX . PDDY . KOUNT . NPR I NT . uC , JC . wY . J Y . DUMMY • 

2CHKY1 .CHKY2.XORIG. YORIG 
COMMON  DELX.ULLY, ID -LX.TDELY.YCOUNT 

C 

RIT  7.  1 00. N. BARL 
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NN*N-  1 

ft  I  T  7  .  1  03. A  I  .  UR2  ♦  I  ft,. 

RIl  /.lOJ.WAUO 

ft  1  T  7. 103* <FORCY(  I ) . 1 *2 .NN) 

-19  0  ftir  7.  1C3.  (  JET  YY(  I)  *  l  uNM) 

DEFINE  INITIAL  STRUCTURAL  CONFIGURATION 

ft  1 T  7.1C2.SIZVEL 

RIT  7  *  103. DEL r IM. TF INAL 

RIT  7 . 1 OO.NF IbRE. ALLRR 

RJT  7 •  l 02 • D «  DPR IM  *  AS 

RIT  7.  I  03 .EC . CEPS Y  «  CEP  SCR  *  CEPSU «  CEP  SO 

ft  I T  7* 1  03.ES ♦ SEPOY 

RIT  7 i 101 . (PAUFAl I ) < 1=1 «NN) 

DO  10  I  *  1  . NN 

10  PALFAC  1  )  =PAUfA( I)*. 01745 

PX(1)«  , 

PY (  1  J  * C. 

00  11  I  =2.  ft 

11  =  1-1 

PX (  I ) =PX(  I  1  )+OARL*COS(PALFA (II)) 

11  »Y(  I  )*PY(  1  I  )+BARL*SlNCPALFA(  II)) 

DO  12  1  *  1  .N 

x  < i ) =PXI I ) 

Y(  I)rPY(  |) 

XOR  | 0 <  I  )  »  * (  1  ) 

12  YORIG<  I  I  * Y (  I  ) 

RIT  7.101 .B.ALERU 
R I T7. 1 01  .PHEE 

RIT  7.101.  GAMMAS. HGT 
R  I  1  7.100.  NPft  I  N  '■ 

WOT  6.231 
WO  T  6.232 
VO  T  0.222.  J 
.VOT  o.23J 

WOT  0.227.(1 .1=1* N) 

WOT  6 . ?2H » ( X ( 1  ) » I  *  1 .N> 

WOT  6.23a 

WOT  6.227. ( I » 1*1 .N) 

WOT  6.228.  (  Y  (  I  ,  «  I  *  )  « '4 ) 

WOT  6.235 

WOT  6.236 

WOT  6.237.U 

WOT  6 . 230 . DECT  I M 

WOT  6.239 

WOT  6.240 

E=HGT#8ARL 

WOT  6.223.E 

E  a  I  728. #GAmma6 

WO  T  6.224.L 

WO  6 . 2  25 

WOT  6.226 

WO  1  6. 227. (  I  • 1 =2 .NN) 

WOT  6.228. (FORCv( I) . 1=2. NN ) 

WOT,  6.229 
WOT  6.230 

WOT  t .c27. ( I .1*2. NN ) 

WOT  6.223. <DtFYY< I ) . I«2.NN) 

V.1 0  T  6  .  .z  3 


J 


j  i  a 


i'l i 


L*#'l~'. 


«0T  6,201 

MOT  *,  203,0. AS 
E«2.*DPRiM 

m  fciiet’£ 

WOT  6  *  206 
WOT  6,207 
WOT  6  ,208 
WOT  6  »  209.ES 
SIGY«ES»SEPSY 
WOT  6,2)0, SI GY 
WOT  6*211 
WOT  6, 209, EC 
SIGY-ECMCEPSY 
WOT  6, 210, SI GY 
WOT  6,212, CEPSY 
WOT  6*21 3  »CEPS>CR 
WOT  6*214  »CEPSU 
WOT  6,215 
WOT  6,216 
WOT  6  *  2 1 7 , A  1 
WOT  6,218«A2 
WOT  6,219. TR1 

»8f  VMi'.m* 

CALL  FTRAP 
CALL  I  NIT 
CALL  DEFORM 
CALL  BET  A 
CC«  1  * 

DO  QO A  I  *  1  ,NN 

pxt i i »x( n 

PY (  I  )  *  Y  <  I  ) 


.j^|  ‘  J|« 

■ 


-.T*  y  <>  >! ;  :i. 

fW  41  ,*  .rf 


'II 


'  i  '1 


■  p<  ■'*! 

train 

t » I  k< ...  1 1 


'  Willi  i  'il  Iff 


004  ALFA  <  I  ) *PALFA (  I ) 

802  IF  < TFINAL-TIME)  800,800,601 
pOl  CALL  LOCATE 
0*  1  • 

call  oeform 
CALL  EXCITE 
CALL  BETA 
GO  TO  802 

100  FORMAT  (  I3.E12.4) 

101  FORMAT  (5E12.41 

102  FORMAT  ( 3E 12*4 ) 

103  FORMAT  (SE12.4) 

200  FORMAT  UH0.10X.19H  SECTION  PROPERTIES) 

201  FORMAT  (11X.19H  **#♦**#  ##«•#***#*,//) 

203  FORMAT  (27H  THE  SECTION  HAS  A  DEPTH  0F.F6.2. 14H INCHES  *N0  M4S. 1*7*4 
1 , 37HSQ  IN  OF  STEEL  PER  INCH  WIOTH  OF  ARCH) 

204  FORMAT  (1H0.49H  THE  STEEL  IS  EQUALLY  DISTRICTED  IN  TWO  LAYEHS  ,F 


16.2.14H  INCHES  APART.,/) 


205  FORMAT  ( 52H  THE  ENTIRE  SECTION  IS  SYMMETRICAL  ABOUT  |Ta  CENTER* » / ) 

206  FORMAT  (1H0.10X.20H  MATERIAL  PROPERTIES) 

207  FORMAT  (11X.20H  #**#**•**  #*####*###  *//) 

208  FORMAT  (1H0.1EH  REINFORCING  STEEL) 

209  FORMAT  (23H  MODULUS  OF  EL AST  I C  1  T Y» , F 1 2 .2 *3HPS 1 ) 

210  FORMAT  l  1  4H  YIELD  STRESS«,F12#2, 3HPS I ) 

211  FORMAT  (1H0.9H  CONCRETE! 
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.ilc  FOR'lAT  (I-+H  V  |  t_LD  -TkAlNa  ,t‘y,5,SH|N/lN) 

13  FOR-'AT  (  1  7h  ORITICAu  STRA  I N*  ,  F9.  5 , 5H  I  N/  J  N  ) 

214  IOHMAT  (  I  7H  ULTIMATE  3  T  H  A  I  Ns  ,  F9. 5  ,  SH 1  N/ I N  ) 

2  1  o  "OHMAI  !  1  HO  «  1  Qx  ♦  26M  PRESSURE  PULSE  PR0PLRT1E3J 

316  OH'm.'  (11X.26H  <»#**  »*i«  sm* ##•#»****«#*#**#■*,// j 

317  (41  P  1  =  ,F9.2  ,  3- iMs  I) 

3  1  ’3  F  OEM  A  I  (AH  P2  =  »  F9  ki  JHPS  I  ) 

219  rOWAT  (1JH  M->L  T I  ML  1  3  .F9.5«7HstC0NDa> 

2*.!  0  f'  URM A  T  (  1  fert  ARRIVAL  T  I  ml  3  3  »  F 9,  s  ,  7  m 3 ECONO a  ) 

?21  FORMAT  tl,*H  <IjL  T  1  ML  2  »  .  F  V.  6  •  7HSEC0N0S  > 

22.?  -'UHMAf  (291 1  Tnt  awlh  13  APPkCX  ImATlU  EY  iUiOH  fclAka,  « /  > 

a2->  1-vRMA  r  (19H  A  RAO..AL  JLR  T  (  Oh  .F7.2.52H  INCMlS  OF  SOIL  la  AS2>OwLO 

i  ro  rejPond  with  the  arch.) 

22*-'  FOl^'Ai  (£6H  ThL  301L  HAS  A  OLNSlTy  0F,r7.2.5H  PCF../J 
22a  )•  yRMA  i1Ja,22H  SOIL  ST RENoTh-Pj l / 1 N , ) 

326  I  I  0  a  t  1  4  H  ****«»******#,/) 

22  ’  •  OH  *R\  [  (  2:1  1  =  ,  13,9!  1  J) 
c2»j  1  OHM  A  r  <  1  IF  1  0,2  ) 

22V  •“OR. MAT  (1HQ.23H  TlELO  DEFORMATION  OF  3PR  I  N&5~  1 N,  > 

221  FORMAT  (1H  «  1  OXi  25H  3  TRUC  T  URAL  CONF  l  CiO’RAT  1  ON*  /  ) 

■  3  FORMAT  (  Im  ,10X,25H  **■!>*»*■»*#*  *«*«•«»***  ) 

S'  F.jpmAT  (  14H  a  C00RO1NATL  ,/) 

234  FC7MA'  (14H  T  COORO  I NA  TL  S  «  /  ) 

-  *  OR  >1 A  r  (  1  HQ  ,  2  r:H  NUMERICAL  INTEGRATION,/) 

2  3c  -  Or  v  a  i  (  1h  ,22h  **#*♦#***  **#*#*#*#*#,/) 

.?_•*’  FORMAT  :6f  Pr  Tlt,c6,4) 

2JS  FORMAT  ( 1SH  T l ML  INltRv-L- ,T  10,6,0H  SLCuNUS) 

23  “J^AT  (1H0.10X.16H  SOIL  PROPER  T  1 F  S,/ ) 

24  0  ;  (llx.lt-H  .►*** 


**  jMltoBtilat >ln 


dimension  and  common  ST atements  on  FJRST  page  of 
COMPUTER  PRINT-OUT  ARE  TO  BE  ADOED  HERE. 


subroutine  beta 

130  FORMAT  <  E  1  0 . 4 , 1  *  ,  E  l  0. 4  > 
132  FORMAT  (E10.4) 


SECONDS *✓  J 

DEUX ( I >  OELY ( I )  OX ( I 


120  FORMAT <//, 13H  T|M£.#E10.».13h 

121  FORMAT  <91M  |  X ( 1 )  V| | > 

*  >  DV( I)  OUX( 1  )  00 Y ( 1 )  > 

122  format  <ij.ix,Eio.4,jx.eio.4.ix,Eio.4,ixEio.4,ix.eio.4.ix,eio.4,ix 

1 .El 0.«. t  X.El 0.41 

fff  ZJSHZ'9**  '  F0WCan>  *»<»>  HOmjtn  DEFORM  BETA  , 

»2*  FORMAT  (91H  *  **••  **••  *•••*•*  •••*•«*  **** 

*****  ******  *mh«  ,/i 

125  FORMAT  (56H**  **********  *******  ********* 

1  *  /  ) 

126  FORMA  T  < / »  38X  * 16HNUMBER  OF  CYCLES) 

127  FORMAT  (  32H  I ‘FVl .FV2.VV1 .VV2.FY.BLASTY  ) 

12S  FORMAT  <  31 H  I .FhI .FH2.VH1 .VH2.FX1 .FX2  ) 

129  FORMAT  <  2 1 H  I .BLASTl .BLAST*  ) 


****** 


acouNT.o. 

DO  5  1*2. NN 

TDELY  (  I  )  ssTDELy  (  I  )  “DEL  Y  <  I  ) 

TDELX (  I  1 «TDELX< I ) -DELX( I ) 

5  CMASSf I ) =<HGT*(BARL**2, J *GAMMAS+ . 0B69*BARL*D » /306 . 4 
45  00  10  1=1 .NN 

FV( I ) «FORCB( I )*S I  K  ( ALFA ( I » ) 

FH ( I ) »F ORCB ( I ) #COS ( ALFA ( I) ) 

VV(  I)  «=VB(  I  >  *COS  (  ALFA  (  J  )  ) 

10  VH(  I)aVB(I  ) *Sl N( ALFA (  I  )) 

DO  20  I»1 .NN 

IF  (  ALFA  C I } |  1.2.3 

1  F  V 1  ( 1  ) «-FV(  I ) 

FV2( I )«FV( 1 ) 

FHI  <  I)  ■  «*FH  (  1  ) 

FH2< I ) «FH( I ) 

VV1  <  I  )  a  VV  <  I  ) 

VV2( I ) a-VV< ( ) 

VH 1(1) «-VH( | ) 

VH2 ( I ) «VH( I > 

GO  TO  20 

2  FV 1(1) >0, 

FV2( I ) *0. 

FHI < I ) «-FH( I ) 

FH2<  I ) »FH(  1  * 

VV1 ( I) «VV( I ) 

VV2( I J  =-VV(  l  ) 

VH1 ( n  30. 

VH2(I)aO, 

GO  TO  20 

3  FV1  (  |  J a  — FV (  I  ) 

FV2(  I  ) aFV (  I  ) 

FHM  I  >«-FH<  1  ) 


,i,  mu  nil-lit.. Ii i<l  'll  nni uni. 


<1  ill  4>lt' 


FH2( I ) *FH< I > 

vv i  ( i )  »w<  n 

VV2  (  I  1  =  -VV<  I  ) 

VH1 < I) =  -VH< I ) 

VH2 (  I  J  =VH<  r  ) 

20  CONTINUE 

BCOUN  T  = BCOUNT  +  1  . 

IF  (  BCOUNT-a.  )  309 , 305 « 303 

305  T|ME»TIME  -  deltim 
DO  308  I  *  1  .N 

WFRT  X  <  l  )«WFRTx( I  >+SJ ZVEL*OELTIM 

XD I FF  X (  I  > ■ Y (  I ) -WFRTX<  I J 

IF  (  XDIFFX(I)  )  306  *  306 «  3C7 

306  T<  I)  «  0 . 

GO  TO  308 

307  T(  I  )*XDIFFX(  I  )/SI ZVEL 
305  CONTINUE 

DELTIM.. 5*DELTJM 
RETURN 

309  DO  30  I=2.NN 
r  i  *  i  - 1 

6  0  FV (  I  ) s-FV2 <  I  I  J-FVl (  I  )-VV2<  I  I  ) - VV 1(1) +  F  Y (  I  )  +BLAST Y (  I  > 

FHtl)  «-FH2  :  !  I  )-FHl  (  [  )  -VH2  <  I  1  J-VHI  {  J  )  — FX  1  (  I  )  +FX2  (  I  >+bl.AST  t  1  l  >*SLAJiT 
12(1) 

PDDX (  I  )  «ODXEND ( I ) 

PDD Y  <  I  )  »DDYEND ( I ) 

DDXEND (  I  )  sFH (  I  ) /CM ASS ( I  ) 

Dl>YEND<  I  )*FV(  I  )/CMASi<  I  ) 

30  CONTINUE 

DO  25  l*2.NN 

AQ  =  ABSF ( ALERS*ODXEND<  I  )  > 

AX  *ABSF ( PDDX (  I  ) -DDXENO  <  I  )  ) 

IF  (  AQ-AX  )  21,22.22 
2  k  IF  <  AX", l  )  22,22,2 6 

22  AQ= A0SF < ALERB*OOYEND ( I ) ) 

Ay.ABSF (PDOYf I )-DDYEND( I ) ) 

IF  (  AQ-AY  )  24,25.25 

24  IF  (  AY-,1  )  25,25,26 

25  CONTINUE 
GO  TO  42 

26  DO  23  1*2, NN 

DELX (  I  )»DX<  I  )*DELTlM+( ,5-ai*DDX( I  > * ( DELT I M*#2 . ) +o*DDXENO <  I >*(DELTI 
1 M*#2. ) 

DXENQ (  1  )*DX(  I  1 +  • 5*DDX  <  I ) *DELT I M+ , 5*D0XEND  { I >*DELT |M 

DELYt  I  )  *D  Y  <  I  )*DELTIM+(  ,5-B)*DDY(  I  >  *  (  DELT  I  M##  =  ,  )  -» U*DDYEND  (  I  )*<UELT| 

I ) 

23  DYEND<  I  >«DY<  I  )+.5*DDY(  I ) «DELT I M+ , 5#DDYEND ( I ) #DELT I M 
DO  400  1=2, NN 

T DEL Y (  I )«TDELY(  I  )+OELY<  I ) 

400  TDELX<  n.TDELXt  I  )+DELX(I> 

41  CALL  DEFORM 
CALL  EXCITE 
DC  401  I =2 ♦ NN 

TOELY (  I ) «T0EL Y  < I ) -DSLY (  I ) 

401  TDELX(  I)*TOELX(  I  )-DELxm 
GO  TO  45 

42  KOUNT  «K0UNT+  1 

00  402  1=2, NN 

TDELY (  I  )  aTDEL Y (  I  ) +DELY (  I ) 


TDELX  <  I  >  -TDELX  <  t  wDELX  (  I  ) 

X( 1 ; «X<  J  )+OELX< I  ) 

402  Y( I ) ■ Y (  I  >+DELY  <  I  ) 

IF  (  NPRJNT-KOUNT  )  43,0*46 

43  WOT  6*  120*  TIME 
WOT  6*121 

WOT  6*124 
DO  SO  l «2  *NN 

50  WOT  6.  122.  1  *X(  I  ,  ,  Y  (  I  )  « TDELX  <  I  )  ,TDELY<  I  )  .DXENDC  I  >  ,£>YEND<  I  )  .UDXtNDl  1 
1  )  .ODYENDl I  ) 

KOUNTrO 
WOT  6*126 
WOT  6*  123 

WOT  6*125 
DO  SI  I ■ 1 «  NN 

51  W0T6.  122 . I ,FORCB<  I ) * V8 < 1)  « QMOM J ( I ) . COUNTS ( I >  *BCOUNT 

46  IF  <  BCOUNT-4,  )  310,44,44 

310  DCLTIM»DELTIM*2. 

44  CONTINUE 

DO  75  1*2, NN 

6U0  PX I  I ) «  X  < I ) 

PYU)iY(I) 

DX<  II -DXEND  (  I  ) 

DDX  ( I >«ODXENO< I ) 

0Y( I) -DYENDl I ) 

DD  Y (  I  ) ■DDYENOI  I J 

DELXI  I  )  =»DX  (  I  )*OELT  I  M4-<  .5~B)*DDX(  I  )  *  <  DELT  I  M*#2,  >+d*POXEND<  I  >*<U£LTI 
1M#*2.  ) 

DELYI  !  )«0Y(  I  )  *OELT  I  M+(  ,5-8>*DDY(  I  >  *M  DELT  IM##2,  )+B*pOVEND(  I  >*<UELTj 

1 M*#2 • » 

DXENO  <  I  )  *DX  <  I  M-.5*UDX<  I  )  *Q£LT  I  M+  .  5*DDXEND  (  I  >*D£LTIM 
DYEND (  I >  »Q  Y (  I )  +  ,5*DDY( I ) *DEL T I M+,S*DDyEND < I ) *DEL  T I M 
TDELY(  I  ) -TDELYI I ) +OELY ( I ) 

TDELX (  I  > a TDELX (  I  >+DELX<  I ) 

73  CONTINUE 
CALL  YIELD 
CALL  CHECK 
RETURN 
END 


Q I  MENS  I  ON  AND  COMMON  STATEMENTS  ON  FIRST  PAGE  OF 
COMPUTER  PRINT-OUT  ARE  TO  BE  ADDED  HERE. 

SUBROUTINE  CHECK 
YCOUNT  =  0* 

DC  13  1*1. NN 

L  C  <  I  )  *  0 

DO  10  I  I »1  ,  NF 1 0RE 

IF  <  TFEPSb  (  Ml)  J  1  1  ,  1  0<  10 

11  IF  (  TFEPSB (1,11)  CEPSu  )  12.12.10 

12  YC  OUN  T b  YCOUN  T+ 1 » 

L  C  (  I  >  a  -  1 

GO  TO  13 
1C  CONTINUE 

13  CONTINUE 

DO  2  3  I  *2,NN 
JC< I ) =0 

DO  20  11*1  .NF IBRE 

IF  (  T FE°3  J (  I , I  I  )  )  21 .20.20 

21  IF  <  TFEPSJ ! 1 . I  I  >+ CEPSU  >  22.22.20 

22  YCOUNT«YCOUNT+l , 

JC( I ) a-1 

GO  TO  23 
20  CONTINUE 

23  CONTINUE 

DO  31  1*1, NN 

YC«0  • 

EPSNEG*0 • 

WWSf  1  .NFI8RE 
IF  (  TFEPSb (  1 . 1  I  )  )  39,38,38 
39  EPSNEG=EP3NEG+1 • 

30  CONTINUE 

DO  30  11*1 .NFIBRE 

IF  (  TFEPSBCI.II)  )  32.30.30 

32  IF  <  TFEPSB! I . I  1 J+CEPSY  )  33.33.31 

33  YC* YC+  l  • 

IF  (  EPSNEG-YC  )  30,34.30 

34  YC  CUNT  *  YCOUNT  + 1 . 

LVI  I ) *  - 1 

30  CONTINUE 

31  CONTINUE 

00  41  1*2, NN 

YC  *  0. 

EPSNEG*0. 

JY ( I ) *0 

D04S  I  I *1 .NFIBRE 
IF  <  TFEPSJ (1,11)  )  49.48,48 
49  EPSNEG*EPSNEG+1  . 

48  CONTINUE 

DO  40  I  I  *1  .NFIBRE 

IF  (  TFEPSJ!  I*  II)  )  42.40.40 

42  IF  (  TFEPSJ! I . II l+CEPSY  )  43.43.41 

43  YC* YC+  1  . 

IF  (  EPSNEG-YC  )  40,88,40 
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80  YC0UNT*YC0UNT+1, 

JY< I )  «-i 

40  CONTINUE 

41  CONTINUE 

44  00  50  I»1 «NN 

IF  (  L  Y (  J )  )  SI « 80*50 
SI  WOT  6 • I  00 • I 
50  CONTINUE 

00  60  l-l.NN 
IF  (  J Y ( I )  )  *1.60.60 
61  WOT  6* 101. 1 
60  CONTINUE 

45  00  70  I ■ l .NN 

IF  <  LC<  I  »  )  T! #70*70 
71  WOT  6.102*1 
70  CONTINUE 

00  80  1-2. NN 

IF  (  JC< II  I  81 *80.80 

81  WOT  6.103.1 
80  CONTINUE 

IF  <  YCOUNT-2.  )  90.81 *91 
91  CALL  SYSERR 

100  FORMAT  <  1 9H  YIELD  HINGE  IN  MR.  13) 

101  FORMAT  <21H  Y IfLO  HINGE  AT  JOINT. 13) 

102  FORMAT  (16H  CRUSHING  IN  BAR, 13) 

103  FORMAT  «  1 6H  CRUSHING  AT  JOINT. 13) 

90  RETURN 

END 


a* .  Maaiim  MOM 


DIMENSION  and  COMMON  STATEMENT j  ON  PJRsT  pace  of 
COMPUTES  PS  I  NT-OUT  A«£  TO  HE  ADDED  HEME* 

SUBROUTINE  DEFORM 
call  FTOAP 

DO  11  I ■ 1  *  NN 

IF  <  ABSF  (Y(  1  +  1  >-Y(  I  M-ABSF  (X<  I+1)-XC  I  )  >*667.96  )  4,j,J 

3  IF  <  Y  <  I ♦ 1 ) « V  (  I  >  )  5.6*7 

5  ALFA ( I )  a-1 .37080 
GO  TO  10 

6  ALF A (  I  )  «0.0 
GO  TO  10 

7  ALFA (  I  ) s  1.37080 
GO  TO  10 

4  AL  f  A  (  J)=ATAN(<Y<I+l)-Y(I))/<X<I*n-X<l))) 

1 0  OUAHL  <  1  ) «<  TOLL  X ( l ♦ 1 ) -TDELX <  1 ) )*COS<  ALFA  I  1  )  )  * ( 7DELY (  I  +  1 I-TOCLY I  I >)* 
1  S I N ( ALFA (  I  )  ) 

1  1  DALF  A  (  I  )  *(  <  TOt£LX<  I-M  )  -TDE!  X<  I)  J  *S  I N  (  ALFA  (  I  >  )-(  tDELY(  J  ♦  1  l-TOELYI  I  I  ) 
1 *COS ( A  LF  A (  I  )  )  ) / (OARL  +  DBARL C  I  )  1 
DO  20  I-t.NN 

20  TEPOARI  I  J«DBARL<  I 1/9ARL 

DC  21  1=2, NN 

DANJ  T  (  I  )=DALFA< I  )-OALFA<  I  —  1  > 

21  CRV~>T  <  I  )  =DANJT  (  I  )  /3ARL 

35  CONTINUE 
CRY j  T ( 1  )  =0 . 

CRVJT ( N) «0. 

DO  36  I ■ 1  «  NN 

36  CRVBARI  I  )«(CSVJT(  I )+CRVJT(  I  +  I j  )/2, 

D°  37  I,2*NN 

37  EPSJTI  I  )  =  ( TEPBAR (  1-1  l+TEPBAR ( I  )  )/2. 

fipre-nfibre 

DO  AO  I«l.NN 
DO  4Q  I  I al  .NFIBRE 
XI  I ■  I  I 

4  0  TFE  PSD  <  1  «  I  I  ) 3  TLP8AR ( I ) - ( D/2. - <  ( ( 2. *X I  1 -1 • ) /2. )*<D/FIBRE> I >#CRVLA«< 

1  I  ) 

DO  41  I  »2  ,  NN 
DO  41  I  I al  .NFIBRE 
X  I  I  «  I  1 

4  1  TFEP3J  <  I  •  1  I  ) =EPSJT(  I )-(D/2.- ( ( <2.*X I I-l . )/2. )#< D/FIBRE) ) ) A CRVJT ( I ) 
DO  42  I ■ l , NN 

TSEP6  T (  I  )  aTEPBAR  <  ! ) +DPR IM« CRVBARf  I J 

42  TStPBU  (  I  )  aTEPBAR  (  I  I  -DPS  1  M -x  CS  VbAR  (  I  ) 

DO  43  I  a 2 . NN 

TSEPJTt  I  )aEPSJT( I  J+DPR1M*CSVJT( 1  ) 

43  TSEPJBl I )«EPSJT< I »-DPRIM#CRVJT( I ) 

CALL  EXTERN 

STRESS-STRAIN  WELAT I ONSHIPS  CONCRETE 
00  290  I a| ,NN 
COUNTB (  I  )  a 0 . 

400  COUNTb  <  1  )  aCODNTBI  I  }  +  1 » 


92 


■  -wii —  «i|—  it—iHfn* ( **<  ■ 


on  on 


00  67  1  1 ■!  tNFIBRE 
Y£PSC*CEPSY-PLSETB<  !  ♦  I II 
IF  I  TFEPSB ( I « 1 1 1  I  59 1 66 *66 
99  IF  (  YEPSC-CEPSCR  I  60 • 60 *350 

60  IF  <  TFEPSBI  I  .  II  l  +  YEPSC  I  62.62,61 

62  IF  (  TFEPSBI I  4  I  I l+CEPSCR  )  64.63.63 
6*  IF  I  TFEPSBI I . Ill+CEPSU  I  66.66.63 

350  IF  <  YEPSC-CEPSU  I  351,351.398 

35 1  IF  <  TFEPSBI I . II l+YEPSC  I  352,3524399 

352  IF  (  TFEPSBI 1 . I  1 I4CEPSU  I  664396,334 
335  IF  I  TFEPSBI 1 . 1 1 l+CEPSCR  I  356.61  461 
356  IF  I  TFEPSBI I . 1 1 i+CEPSU  I  66*61.61 

61  FS  IGBI  1  ,  I  I  !i.EC*<  TFEPSBI  I  4  I  I  l-Pl.SETBI  I  «  I  III 
GO  TO  67 

63  FS IGBI  l  ,  I  , )«-EC*CEPSY 
GO  TO  67 

65  FS IGBI  1  , II) « I-EC*CEPSV >* I  I »- I -TFEPSB I  I . II > -CEPSCR )/ I CEPSO-CEPPCR > I 
GO  TO  67 

66  FS IGB I  I  ,  I  1  I »0, 

GO  TO  67 

334  FS IGBI  I  .  II  I ■ I -EC*CEPSy ; *  I  i .-(-TFEPSB! I « I  I  I -CEPSCR I / I CEPSO-CEPPCR I  I 
GO  TO  67 

356  FS IGB.  I  ,  I  I  ) bEC*! TFEPSBI  i , I  I I-PLSETBI  1,111) 

DUMMY  *  I «EC*CEPSY)*< 1 .-( -TFEPSBI I , I  I  I -CEPSCR I / l CEPSO-CEPSCR I  I 
IF  (  FSIGBI I , I ! 1-DUMMY  )  3S7.67.67 

357  FSIGBII.il  I  ..DUMMY 

67  CONTINUC 

CONCRETE  STRESS-FORCE  RELATIONS 
CFOPCB  I  I  ) =0, 

DO  200  I  1*1  . NF I BRE 

200  CFORCBI  I  )  oCF ORCB I  I )  +  ( I D/F I BRE )*FS I GB ( I .III) 

TOP  REINFORCEMENT  IN  BARS 
SCRATO«ES/EC 

3EPYPT1  J  )»SEPSY+SPSTBT( I ) 

SEPYPC I  1  )«-SEPSY+SPSTBT( I ) 

ESEPBTI I ) «TSEPBT  < I I-SPSTBTI I ) 

225  IF  I  SEPYPCI I )  J  226.232.233 

226  IF  I  TSEPBTI II  )  227 .230 .23 l 

227  IF  {  TSEP8T I  I I-SEPYPCI I )  >  236.238.226 
226  IF  <  SEPYPTi:)  )  229.240.240 

229  IF  I  TSEPBTI  I l-SCPYPTl I  I  )  240.239.239 

230  IF  I  SEPYPT I  I  I  I  239.239.240 

231  IF  I  SEPYPT ( I )  )  339.239.233 

232  IF  I  TSEPFTdl  )  238, 238.234 

233  IF  (  TSEPBTI I ) -SEPYPT ( I )  )  240,239*239 

234  IF  I  TSEPBTI I I-SEPYPTI I )  J  240*239*239 
233  IF  I  TSEPBTI I )  )  239*239*237 

236  IF  I  TSEPBT ( I I-SEPYPTI I )  )  240*239*239 

237  IF  I  TSEP8T  t ] )— SEPYPCI I  I  I  236*238*236 

238  SFORBTI  I  ) »-ES »SEPSY*I SCRATO— 1 • J *AS/ 1 2#*SCRAT0 ) 

GO  TO  241 

239  SFORBT I  I  )  «ES#SEPSY*AS/2 . 

GO  TO  241 

240  IF  I  ESEPBTI 1 )  )  242,242,243 

242  SFORBTI  I  ; «ES*ESEPBT I  1 ) * I SCRATO-i • I *AS/ 1 2,*SCRATO) 

GO  TO  241 

243  SFORBTI  I  )  *ES#£SEP8T I J  )*AS/2, 
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n  r 


241  CONTINUE 

C  BOTTOM  REINFORCEMENT  IN  BARS 

S£PYPT< 1 >*SEPSY+SPSTB3( 1 ) 

SEPYPC (  I ) ■-SEPSY f$PST0B  < I } 

ESEP88 ( 1 ) "TSEPBBf 1 ) -SPSTBBI 1 > 

325  IF  (  SEPYRC(I)  )  326.j32.335 

326  IF  (  TSEFB9U)  I  327*330.331 

327  IF  (  TSEPB0!  I  I-SEPYPCl  |  }  |  336*338*326 
320  IF  (  SEPYPT(I)  )  329*340*340 

329  IF  {  T5EPBB1  I  )-SEPYPT<  I  1  j  340*339*339 

330  IF  (  SEPYP T  <  !  )  >  339*339*340 

331  IF  (  SEPYPT(I)  )  339.339,333 

332  IF  (  TSEPBB(l)  )  338.338.334 

333  IF  (  TSEPOB!  I  I-SEPYPTi  I  J  )  340*339*339 

334  IF  <  TSEPSBl  I  I-SEPYPTt  I  )  )  340*339.339 

335  IF  (  TSEPBB(I)  )  338.336.337 

336  IF  (  TSEP8BI  I  >-SEPYPT(  I  |  )  340.339.339 

037  IF  (  TSEPBBI  I  J-SEPYPCI  1  )  |  338*338.336 

33B  SF QRBB  <  I  ) »— E J*SEPSY* I  SC RAT 0-  l  # ) *AS/ < 2 . *SCHATO ) 

GO  TO  341 

339  SrORUa ( I) »ES*SEPSY*AS/2. 

GO  10  341 

34j  IF  (  CSEPBB(l)  1  342 , 342 .343 

342  SFORBBI I  >  =ES*ESCPBB< I  ) *<  SCRATQ-l t )*AS/(2.*SCWAT0) 
GO  TO  341 

343  SFORBB(I)  *ES*£SEPBB< I )*A5/2. 

341  CONTINUE 

F OWCB  <  I  ) *CF QRC0  < i ) +SFORBT (  I I+SFORBB! I ) 

IF  <  r^EPSB  I  I  .  I ) *TFtPSB< 1 .NFIBREI J  801*801*290 
001  FORUFB!  I  >«FORCem-EFO«CB<  n 
oqkabsfifordfbi i i j 

IF  <  OQ-AL.ERR)  290*290.293 

293  IF  <  FORDr O (  1  > )  294,290*295 

294  DIKLPBI I )=FORDFB< 1 ) / < EC* ( 0+ ! SCRAT 0- I . > *AS ) > 

GO  TO  296 

295  OIFEPBI I )»FORDFB( l)/(ES*AS) 

296  DO  297  I  1*1 ,NP1BRE 

297  TerFPSB(  I  .  I  I  )  "TFEPSB  (  I  .  I  I  I-DIFEPB!  I  ) 

TSC°BT< I >»TSEPBT< 1 ) -O I FEPB ( I > 

TSEPBBf  I ) *T  SEPOB ( I J-OIFEPBI I > 

GO  TO  400 
290  CONTINUE 


00  291  1*2*  NN 

COUNTBC I ) *0  » 

450  CO'JNTB  (  2  )  *C0UNTB  (  I  )  + 1  • 

DO  87  I  1*1  .NFIB9E 
YEPSC«CEPSY-PLSETJ<  I « I  I  > 

IF  {  TFEPSJ! 1,11)  )  79,86.86 

79  Ip  (  YEPSOCEPSCR  )  80,80.370 

80  IF  (  TFEPSJ!  I.IIM-YEPSC  )  82*82*81 
B2  IF  (  TFEPSJ  (  I  .  I  IH-CEPSCR  >  84,83,83 
84  IF  (  TFEPSJ (1,11 I+CEPSU  )  86*86*65 

370  IF  (  YEPSC-CEPSU  )  371,371,378 

371  IF  (  TFEPSJ! I , »I)4YEPSC  )  372,372,375 

372  IF  (  TFEPSJ! 1 , I l I+CEPSU  I  06,374,374 

373  IF  I  TFEPSJ!  I  ,  m+CEPSCR  )  376,81  ,81 
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378  I F  (  TFE.PSJ  (  I  «  i  I  1*CEP$U  1  86*81*61 
81  FS  I  G  J  (1*11) -EC4ITFEPSJ! J  .  I  1  1-MLSETJ!  I  i  1  1  »  ) 

GO  TO  67 

83  FSIGJ<  I  .  lI>«-eC»CEI»4Y 
GO  •f,0  87 

85  FSIG--!  I  *  III  - <-EC*CIP»Y 1*(  1 .-!-rFCPSJ!  1*11 1 -CEPSC& ) / < CEPSO-CEPSCR 1 ) 
GO  TO  87 

•  6  FSIGJf  1*11  1  **Q  • 

GO  TO  87 

374  FS IGJ (  I  ,  11 ! ■ < ~EC*CEPSY 1*(  1 (-TFEPSJ!  1*11  I -CEPSCR 1 / < CEPSO-CEPSCR > ) 
GO  To  87 

376  FS1GJ <  1  *  ID «EC*< TFEPSjt i . I  1 ) -PLSETJ! 1*1111 

DUMMY  * ( -EC«CEPSY ) * <  1 . - ( -TFEPSJ  < I « II  1 -CtPSCR ) / <  CEPSO-CEPSCR I  1 
IF  (  FS!GJ<  1  *  I  I 1-OJMMy  1  377*87*87 

377  FSIGJI  I  *  I  I >»OJMMY 
87  CONTINUE 

CF  OWC  J (  I  )«0. 

C  MOM  J (  I  ) «0. 

00  201  11*1 fNFlBBE 
X  I  I  *  I  I 

CFORCJ ( I ) -CFORCJ! I  1 ♦ C ( 0/F I BRE 1 *FS I G J <  1  •  I  I  1  I 
201  OCMOMJ  (  1  )  «*C  MOM  J  (  I  )  F  S  IGJ  <  1*111*1  1  0/F  I  8  RE  1  *  ! 0/2.  -  <  I  «  *•**  1 1  - 1  4 
1D/F1BRE1 1 ) ) 

TOP  STeEL  AT  JOINTS 
SEPYPT!  I  I-SEPSY+SPSTJTI 1 } 

SEPYPC  <  I  I --SEPSY+SP3TJT ( 1 ) 

ESEPjT  (  I  J-TSLPJT!  I  I  -SPSTjT  (  I  ) 


425 

IF 

( 

SEPYPC ( I ) 

I  426*432*435 

426 

IF 

( 

TSEPJT! I J 

)  427*430*431 

427 

IF 

< 

TSEPJT ( I ) 

-SEPYPC! 11  *  438*430.428 

428 

IF 

( 

SEPYPT  (  t  ) 

}  429*440*440 

429 

IF 

( 

TSEPJT  !  1 ) 

-SEPYPT!  |1  )  440*43*'  439 

430 

IF 

< 

SEPYPT  (  I  ) 

)  439*439*440 

431 

IF 

( 

SEPYPT  !  I  1 

)  439*439*433 

4  32 

IF 

( 

T SEPJT  (  I  1 

)  438*438*434 

433 

IF 

( 

TSEPJT!  1  1 

-SEPYPT! 1)  >  440*439*4 

434 

IF 

( 

TSEPJT ( I } 

-SFOYPT(I)  |  440*439*4.7, 

435 

IF 

( 

TSEPJT! I | 

)  438*438*437 

436 

IF 

( 

TSEPJT! 1 ) 

-SEPYPT! I)  )  440*439*439 

437 

IF 

< 

TSEPJT! I ) 

-SEPYPC! I)  1  438s 438*436 

4  38  3F  OR  JT  (  I  )  «-Ei»*S£PSY*  <  SC  PA  T  0-  1  *  }«AS/I  SCRAT  *2*1 
GO  TO  441 

439  5F0RJT1 n «ES*SEPSY#AS/2. 

GO  TO  441 

440  IF  (  CSEPJT(I)  )  442.442.443 

442  SF OR J  T (  J )«CS*ESEPJT(  I ,*< SCRATO-l • I *AS/ < 2.*SC«AT0 1 
GO  TO  441 

443  SFORjTl | |-ES*ES£PJT< i )*AS/2. 

441  CONTINUE 

C  BOTTOM  STEEL  AT  JOINTS 

SEPYPT!  I  I-SEPSY+SPSTJS! I  1 
SEPYPC (  I  1--SEPSY+SPST JB! I j 
ESEPUB ( 1 l«TSEPje( I 1-SPSTUB( I > 

525  IF  (  SEPYPC <  I )  I  526*532.535 

526  IF  <  TSEPJdll)  1  527*530*531 

527  IF  (  TSEPJSC  1  I -SEPYPC!  1 )  J  538.636*528 

528  IF  (  '.EPYPT  (  1  )  1  529*540.540 


55 


529  IF  (  T  5E f-'J-J  (  1  )  -SEP  YP7  (  I  )  )  5*0  *  *39 ♦  339 

530  IF  (  SEPYWT  <  >  )  539  •  539  .  5*0 

531  IF  (  L-r  n  yp  T  (  1  >  )  539.539.533 

532  IF  (  Titojb'I)  )  536*538.53* 

333  IF  (  r  jtPja l  I ) -SEP YP  T ( I )  )  340*339.539 

534  IF  (  TSEPJFM  I l-SEPYPT I  I >  )  540*039*539 

535  IF  (  TSEPJtM  I  )  5  536*538.537 

=  36  IF  <  T.'tPJt:;  (  I  ) -bliPYPT  (  I)  )  540*539.839 

537  IF  (  TSGPj.ti  1  )-SEPYPEt  I  >  >  530.538*536 

530  3FORJ0I  I  ) =-Ea*5EPo Y* ( 5CPAT0-1 , > *  AS/ < 2 « *5C9A TO } 


GO  r:  39) 

j39  UkJU <  I  ) b*iEP5Y*Ao/2 • 

GO  TO  3*  1 

"40  (  ■- i'-PJb  (  I  )  I  542  *  54  2  *  343 

-42  ijFLWJdll  l*L.3t  SEP  Jp  (  1  )  *  (  jCkATO-  U)*Aj/(i.*iiOKATw) 


GO  TO  541 

54  3  _F09Jt3<  I  )*L:S*tiEPj-J(  I  )  •’■A  3/2  • 


CONTINUE 


I  OPT  „■  (  I  )  iCFO^Cj  (  1  )  +SFORJT  (  1  )-*SFORJB<  1  > 
j  i..  _ j  (  |  j  =c  wON-i  (  I  ) +0PRIM* ( SFQRJB I  I  ) -SF09JT (  I  )  ) 

I F  (  TFLP3UI  I  .1  )*TF£PS3(  I  .NFJBPE))  600 » 600  * *9 1 
*•  ,<Z  FUPOFJI  I  )*FG9CJ<  I  I-t.F  Ja'CJI  I  I 
j:)eAi$r  (FOWLiF  j  (  I  }  ) 

IF  (  OU-AuERR)  291*291.451 
•  -1  IF  (  F CRDF  J (  I  )  J  432.291.453 

4  32  10  I  F-.P  J  (  I  )  *FQ«UFJ  (  I  )  /  (  EC  *  <  0+ <  SCHATQ- 1  .  )  *AS)  ) 

GO  TO  454 

4  53  D  I  F  [■;  p  J  (  1  J.FOBDFUI  I  )/  (£3*  AS) 


4  04 
4  'a 


Do  455  I  I 
T 


:  I 


(  I  .  I  I  5  -  T 


SJ<  I . 1  I  I-DIFEPJI 1  ) 


T  SEU  J  T (  j  ) rTSLPjTI  I  l-DIFEPDl  t ) 
T5E°J3<  I ) a T SEP JB (  I  I -DIFEPUI  I  > 

GO  To  450 


291  COT  "  !  .vC 

00  292  1*1. NN 
I  I  <•'  I  F  1 

V3  <  I  )  * (  l */ I 0A9L  fDBARE (  1  )  )  )*IBM0MJ<  I  I  ) -8M0M J ( I  >  > 


o  ■< :  :  on  r  i  n  j- 


9ETUBN 

END 


MLjuumsjm  img 


dimension  and  common  statements  on  first  PAGE  of 
computer  print-out  are  to  ae  added  mere. 


SUBROUTINE  excite 

c  calculate  pressure  amplitude  at  time  t*t 

DO  600  I*2«NN 

IF  (  T(I)-TR»  )  30S.303.306 

305  BLAST  * < T (  I) /TR1 I *A 1 
GO  TO  314 

306  I  f-  (  T  (  I  )  -  T  AR2  )  307.307.310 

307  IF  (  TR2-TAR2  )  308.309.308 

308  BLAST ■ A  1 
GO  To  314 

309  BLAST* A  1 ♦ (  < A2-A1 )/( TAR2*TR1  >)*IT<  I  J-TRI) 

GO  TO  314 

310  IF  (  T  <  I  ) — T R2  )  311,313.313 
31t  IF  (  TR2-TAR2  >  312.313,312 

312  BL  AS  T* A 1> ( ( T ( I ) -TAR2) / ( TR2-TAR2 ) ) *  < A2-A1 ) 

GO  TO  314 

313  BL  AST  *Ae 

314  CONTINUE 
I  I  *  !  -  1 

CHECK  TO  DETERMINE  the  PRESENCE  OF  SPRING  1  ANQ/OH  SPRING  2 

VV«0. 

Z2*0. 

I J* 1+  1 

IF  (  X (  I  )  — x (  I  I  )  j  406,407.408 

406  ZX«-1. 

GO  TO  409 

407  X<  I >  *X  <  I )*.01 
ZX*I  . 

ZZ*  1  * 

GO  TO  409 

408  ZX*1. 

409  CHKY  1  (  1  )  *Y<  I  )  -  <  (  Y<  I  )-Y  (  I  I  )  |oK(  1  I  I  1  >«ZX 

IF  (  ZZ  )  411,411,410 

410  XIII.XII >-.01 
ZZ*  o. 

411  IF  (  X  (  I J ) - x (  I  )  )  4  1  2 , 4  I  3  1 4 1 1 , 

412  ZXS-U 

GO  TO  415 

413  X<  I J ) *  X (  1  J)+.01 

zx«; . 

ZZ  *  1  • 

GO  TO  415 

414  ZX*1. 

415  CONTINUE 

CHKY2I  I)*Y<  I  )  +  (  <  Y(  I J )-Y (  I  >  )/<  XI  I J>-X (  I )  ) ) *2X 
IF  <  ZZ  )  417.417.4J6 

416  X  <  I J  J  *  X (  I JI-.01 
ZZ*0. 

417  IF  (  y (  I  ) -CHK y 1 (  I >  )  400,400,401 

40C  FXl  (  I  >  *0. 

BL ASTI  <  I  ) =C . 


f 


^  "'f  -t  t  # 
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Y  Y  a  1  , 

*01  OEFX2 (  I  ) --T0LLXI  1  > 

IF  <  Au)SF<QEFx2<  1)  >-.0001  )  434.402*402 

*34  0EFX2 ( I ) *0. 

402  IF  <  Y<  I  )-CHKY2< 1 )  )  403.403*404 

403  FX2  <  I  )  «0. 

BLAST  2 (  1)*0. 

VV«1  . 

«04  DEFX1 ( 1 >«T0FLX( 1 ) 

*05  DEFY!  1  )  *-riOELY(  I  ) 

IF  (  A0SF ( UtFX 1 ( 1 ) ) - • 0C0 1  )  430.431,431 

430  DEFX l (11*0, 

431  IF  l  ABSF(DEFy< I )) -.0001  )  432.433.433 

432  DEFY!  I  ) *0. 

*33  IF  (  YY  )  360,365.366 

365  IF  (  DEFXl(I)  1  212.213.213 

212  FX1 ( I ) »(FOWCY( 1 ) /D£FYY< ! ) ) *D£FX1 ( J ) 

00  TO  224 

213  FX 1  ( I  )  *0  • 

22*  auij  r  1  (  I  )  =wflT  1  Q*SL  AST  *  <  ABSFI  •  5*fciARL*6  I  N  (  ALF  A  <  |  1  )  > ,) +nr;>,aF  <  ,a4o**«L*»J 
INI  AuFAl  1  )  )  )  ) 

366  CONT | NUt 

IF  <  VV  )  465.465,466 
465  IF  (  DEFxctn  )  318.31W.319 
3  I  Fx2<  I  )  s(F0HCy(  I  )/DEFYY(  1  )  )*DEFX2<  I  ) 

G--  TO  324 
11  W  F  y  2  (  I  )  s  o  * 

■  4  6LA.J  T2  (  I  )  s-WAT  I0*0LA6T*  I  AB6F  <  .5*e»ARL*S  1  N  <  ALF  A  <  1  |  »  )  )+AuSF(  *&*u»Rl*S 

1  ! M  ALF A (  I  >  )  )  ) 

*66  .oNT  IN^t 

I  F  (  j L .  F  Y  (  1  )  )  4  J  8  »  *  1  9 , 4  1  9 

*10  FYI  I  >  =  (F0RCY1  I  )/DLFYY(  l  M»OEFY(  I  ) 

GO  TO  424 
•19  FY( |)»0, 

*  ,  ulujTV  (  I  )  --JL.AjT*l  AbSF  (  .5*oAnu*LUS(  ALF  A{  111)  >+Ao;»F  (  .6-*fc)ARL*CWa  (  ALF 

1  A  I  I  )  )  )  ) 

066  CONTINUE 

KX  1  (  I  )*FX1  <  1  )  *  (  AB3F  I  «5*t!AWL*S  1NT  ALFA  (  111)  H-AoaFC  .  »*uARL#&  I  N<  ALF  A  (  1 
1  >  )  )  ) 

Fx2(  I  )«FX2<  I  )*(ABSF  l  •3*UARL**6!NC  ALFA  (  11))  H-AUSF  (  .»+UAWL#»lN«  AwF*  (  I 
1  )  )  )  ) 

F  Y  (  I  )  »F  Y  (  I  )  *1  ABSFI  .5*BARL*CG5<  ALFA!  11))  )+A;*SF  <  •  5*  o  A  *L  *  COS  (  ALFA  <  J  )  ) 

1  )  ) 

'.-00  CONTINUE 
•E  TURN 
;.ND 


M 
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o  n 


W<It|  1U|  ll»4nk 


i.  liiottlii  UffM| 


DIMENSION  ANO  COMMON  STATEMENTS  ON  FIRST  PAGE  OF 

computer  print-out  arc  to  ae  added  here. 

subroutine  extern 
DO  67  1  ■ 1 «NN 
EYEPSC»CEPSY-EPSETB( I ) 

IF  <  TEP2A9 ( I |  )  39(66,66 
39  IF  (  EYEPSC-CEP6CR  )  60 « 60 « 330 

63  IF  1  TEPBARI I J4CYEPSC  )  62,62,61 

62  IF  (  TEPBARI I )4C£PSCR  )  66,63*63 

64  IF  (  TEPBARI I I4CEPSU  )  66,66,66 

350  IF  (  EYEPSC-CEPSU  J  351,331,338 

351  IF  (  TEPBARI I )6CTCPSC  )  332,352,335 
332  IF  (  TEPBARI I >«C«PSU  )  66,394,354 

353  IF  I  TEPBARI I ) +CEPSCR  )  336,61,61 
J5B  IF  <  TEPBARI I l+CEPSU  I  66,61,61 

61  EFORCBI I  I -EC* I TEPBARI l J-EPSETBI I >  >*D 
GO  TO  67 

63  EFORCBI  I  1 ■~EC#CEPSY*D 
GO  TO  67 

65  G0°T0BF  7  * ’ ( 'tC#CCPiY 1  * ‘ 1 • - < -TEPBARI  I ) -CEPSCR ) / ( CfPSO~CEPSCR  » ) *0 

66  EFORCB  I  I  ) =0 • 

GO  TO  67 

354  EFORCB I  I  )  > I -EC*CEPSY ) *  I l . - I -TEPBAR ,  I > -CEPSCR ) / <  C£PS0-CIPSCR ) )  *0 
GO  TO  67 

356  EFORCBI  I  )  «EC* ( TEPBAR ( (  l-EPSETBl I  1 1*0 

DUMMY ■ I -EC >CEPSY ) *  I  I , - 1 -TEPBAR I  I ) -CEPSCR ) / ( OEPSO-CEPSCR  J  j  *0 

IF  I  EFORCBI 1 ) -DUMMY  )  337,67,67 

357  EFORCBI  I  ) -DUMMY 

67  CONTINUE 

FORCE-STRAIN  RELATIONS  FOR  CONCRETE  AT  JOINTS 
DO  76  I«*.NN 

76  rSPJTl  i  j  «,5*TEPBAR(  1-1 )+.5*TEPBARl I ) 

DO  B7  I *2,NN 
E Y EPS C» CEPS Y-EPSETJ ( I ) 

IF  I  TEPJTII)  >  79 , B6, 66 

79  IF  i  EYCP3C— CEPSCR  )  80,80,370 

80  IF  i  TEPJT I  I  ) 4CYEPSC  J  B2 ,02,81 

82  IF  I  TEPJT I  I ) -f  CEPSCR  >  84,83,83 
84  IF  I  TEPJT! I J+CEPSU  >  86,86,85 

370  IF  I  EYEPSC-CEPSU  )  371,371,378 

371  IF  (  TEPJT!  n+EYEPSC  )  372,372,375 

372  IF  I  TEPJT! 1 l+CEPSU  J  86,374,374 
375  IF  I  TEPjT (£ J+CEPSCR  }  376*81,81 
J78  IF  I  TEPJT I  I ) +CEPSU  )  86,81,81 

81  EFORCJI  I  1 «EC*<  TEPJT I  I >-EPSETj< I ) )*D 
GO  TO  87 

83  EFORCJI  I  1  a-EC*CEPS Y#0 
GO  TO  87 

85  EFORC J  I  I  >  a  I -EC»CEPSY )* I i «-  I -TEPJT 1  I ) -CEPSCR) / I CLPSO-CEPSCR ) >*U 
GO  TO  87  mymit 

36  EFORCJ  I  I  ) =0, 

GO  TO  87 


n  o 


374  EFORC J  (  I)  «  <  -£C*CEPSY  MM  1  .  -  <  -TEPJT  (  1  >  -CEPSCR  )  /  <  CEPSO-CEPSCR  )  >  *u 
GO  TO  87 

376  EFORCJ<  t  )=EC*!TEPJT! I )-EPSETJ< I > )*D 

DUMMY- ( -£C*CEPSY>*<  l.-l-TEPJT!  I > -CEPSCR ) / i CEPSO-CEPSCR ) 1*0 
IF  <  EFORC-i  !  I  ) -DUMMY  )  377*,  87 »  07 

377  EFORCJ!  I  I-DUMMY 
87  CONTINUE 

C 

C  FORCE-STRAIN  RELATIONS  FOR  STEEL  IN  BARS 

SCRATO-ES/EC 
DO  222  I >1 .MN 
SEPYPT (  I  )  -SEPSY4ESPST0 ( 1 ) 

SEPYPC I  I  )  s-SEPSY+ESPSTB t I ) 

EESEPB ( I ) -TEPBAR  < I ) -ESPSTB  t  I  ) 

IF  (  SEPYPC < 1 )  )  226*232.239 

226  IF  I  TEPBAR ( I )  )  227.230*231 

227  IF  (  TEPBAR! I) -SEPYPC! I)  )  230.238.228 

228  IF  1  SEPYPTII)  )  229*240.240 

229  IF  (  TEPBAR!  I) -SEPYPT  !U  )  240.239.239 

230  IF  (  SEPYPT! I)  )  239.239.240 

231  IF  (  SEPYPTII)  )  239.239.232 

232  IF  (  TEPBAR! I)  )  238.230.234 

233  IF  (  TEPBAR! I ) -SEPYPT! I )  )  240*239.239 

234  IF  (  TEPBAR! I )-SEPYPT| I )  )  240*239.239 

235  IF  !  TEPBAR! I)  )  238*238.237 

236  IF  (  TEPBAR! I J-SEPYPT! I )  )  240.239.239 

237  IF  !  TEPBAR! I )-SEPYPC< I )  )  238.238.236 

238  EFORSB!  I  )  =-Ei*SEPSY* ! SCR A TO-1 . )*AS/SCRATO 
GO  TO  241 

239  EFORSB!  I  > *ES# 3EP5Y*AS 
GO  TO  241 

240  IF  (  EESEPB!!)  >  242.242.243 

242  EFORSB! 1 ) «ES*EESEPB ( I ) * ! SCRATO-1 • ) »AS/SCRATO 
GO  TO  241 

243  EFORSB!  I)*ES*EESEPB!  l)*AS 

241  CONTINUE 
222  CONTINUE 

FORCE-STRAIN  RELATIONS  FOfl  STEEL  AT  JOINTS 

DO  262  !*2.NN 

SEPYPT!  I  )*SEPSY+ESPSTJ! I ) 

SEPYPC I  I ) ■ -SEPSY+ESPST J ! I ) 

EESEPJ  !  1  )«TEPJT!  D-ESPSTJ!  I) 

IF  !  SEPYPC! I)  )  326.332.335 

326  IF  !  TEPJT ! I )  )  327,330*331 

327  IF  !  TEPJT! I J-SEPYPC! I )  )  330.338,328 

328  IF  (  SEPYPTII)  )  329,340.340 

329  IF  (  TEPJT! I ) -SEPYPT! I )  )  340.339.339 

330  IF  !  SEPYPT ( 1 )  )  339.339.340 

331  IF  !  SEPYPTII)  )  339.339.333 

332  IF  <  TEPJT! I)  )  338,338.334 

333  IF  (  TEPJT!  I  J-SEPYPT! I  )  )  340.339.339 

334  IF  (  TEPJT! I ) -SEPYPT < I)  )  340,339.339 

335  IF  !  TEPJT! 1)  )  338.338.337 

336  IF  I  TEPJT(  D-SEPYPT!  I  )  )  340.339,339 

337  IF  !  TF.Puf!  1  J-SEPYPC!  I  )  >  33B.338.336 

338  EFORSJ  '  I  )  «-ES#SEPSY# ! SCRATO-i . >*AS/SCRATO 
GO  TO  341 

339  EFORSJ  <  I  ) -f.3*SEPSY* AS 


60 


- lU*l  JW4ki«4iU(»(lWajUi.uJUjiuL,  uUUil 


. „;!i:ii;.ij,  ;■ .  - 


GO  TO  3*1 

3*0  IF  (  EESEP-M  I  I  »  3*2*  3*2 « 34 3  ?  s,  (  | 

344.*  dFOBS>'<  I  )*tS>*ttSEPJ(  «  |*<SO-ATO-1  .  >«/»S/SCMHTO 

GO  TO  341 

3*3  EFORSJl  l)rfcS«EESePJ<  1»«AS> 

3*1  CONTINUE 
262  CONTINUE 

SUMMATION  OF  THRUSTS  IN  BARS  A NO  JOINTS 

DO  300  1 • t • WN 

300  iforcsi  i  > -crones m*EFO«5B;  1 » 

DO  301  I *2  «NN 

301  KFORCJl  I >»EFORCJ( 1 )+EFORSJ<  I  1 

RETURN 

END 


DIMENSION  AND  COMMON  STATEMENTS  ON  FIRST  OF 

COMPUTER  PRINT-OUT  ARE  TO  EJE  ADDED  MERE* 


SUBROUTINE  I  NIT 
DO  76  1*1 «NN 
tepbar<  n*o. 

PBARL < I ) "BARL 
SPSTBT( I ) =0 • 
jPSTBB( I  ) *0  « 
TSCP9T( l ) “0 • 

T  SK.P8B  <  11*0. 
EPSETL1  (  I  )  =0. 

LaPj TB (11=0. 

UO  76  I  I  *  1  .NK I dRE 
T.-EPSB  (  I  .  n  )  *0* 

UL  SL  rB (  I  . I  1) *0« 

7 to  CONTINUE 

JO  78  I  =  1  .N 
TSFPJT ( 11=0. 

tsepjb < I ) =  0. 

iPSTjK  I  1=0. 

jP S  T  JB <11*0. 

3MM3lli:8: 

DO  73  11*1  « NF I 3HE 
Tr'LPOJ  (1.111*0. 

PL SET j (  1 . I  I  1 *0. 

78  CONTINUE 
CC»0. 

T  I  ME  =  0 . 

<OUNT=NPRI NT 

0*0  . 

DO  6  I*2.NN 
DELX< I ) *0. 

OELV( 1 1*0. 

TULLX ( 11*0. 

TDELYC  I >  *C • 

EPSET  J  < I )»0. 
ESPST J ( 11=0. 
ODXEMiM  11*0. 

JO YEN J (11*0. 
u/ENDC  1  1  *0. 

JYCNO ( 11*0. 

DDOX< I) *0. 

PDDY ( 11*0. 

3LA3T1 ( J 1*0. 
BLASTS (11=0. 
BLASTYI I ) *0 • 

*X1 ( I >  *0. 

FX2(  I  1  *0. 

*  y  <  n  *  o « 

DDX ( 11*0. 

UO  Y  (  I  1  *  0  . 

DX< 11*0. 

6  DY( 1 1 *0. 


DIMENSION  ANO  COMMON  STATEMENT..  ON  FIRST  PAGE  OF 
COMPUTER  PRINT-OUT  ARE  TO  BE  ADDED  MERE. 

SUBROUTINE  LOCATE 
IF  (  Q)  30S«309«310 
J«1 

DO  300  1*2. N 

IF  <  Y(J)-Y(I»>  *99.299.300 
J*  I 

WAVE  FRONT  AT  J 
CONTINUE 
DO  304  I  *1  »N 
WFRTX  <  I ) »V ( J ) 

LOCATION  OF  WAVE  FRONT  AT  TIMET+OCLTA  T 
T I ME*T I ME+DELT IM 
DO  308  I *1 «N 

WFRTX ( I ) »WFRTX( 1  I  —  S I ZVEL*D£LT1M 
XD I FFX  <  I) «Y(  I  J -WFRTX (  I  I 
IF  (  XDIFFX(I)  )  306.306.307 
T  <  1  )  *0* 

GO  TO  308 

T  (  I  )  *XOIFFX(  I  J/S1ZVEL' 

CONTINUE 

RETURN 

END 


..  w!  .tmi  t  tpf»  ■ 


•  II  f  I  Irj  I 


DIMENSION  AND  COMMON  STATEMENTS  ON'  FIRST  PACE  OF 
COMPUTER  PRINT-OUT  ARE  TO  BE  ADDED  HERE* 


SUBROUTINE  locate 
IF  (  0  »  305 « 305. 3 1 0 

305  J*1  i1*' 1  ■ 

DO  300  I«2«N 

IF  (  V  I J ) - Y (  I >  )  299 • 299. 300 

299  J«  I 

300  CONTINUE 
PHEE*PHEE*.017A5 
DO  30 1  I « I «N 

WF  RTx  I  1 >  =-l 0000000.0  +  Y ( 12*1 • /< S I N C PHE6 ) /COS ( MwtC > » 

301  XDIFFX! 2  >«rAQ3F  (XI  1  J  -WFRTXI  I  >  ) 

I  I  »  I 

DO  303  I*2.J 

IF  (  XD I FFX ( I ) -XDIFFX { II)  >  302.302.303 

302  I  I «I 

303  CONTINUE 
0=1  • 

WAVE  FRONT  AT  MASS  II  AT  T I M£  T»0 
DO  304  I  *  I  . N 

304  WFRTXI  I  ) bWFRTXI  I  )+XOIFFX(  I  I ) 

LOCATION  OF  WAVE  FRONT  AT  TIME  T  «Q+DEC  TA  T 
310  TIMEsTIME  +UtL  T I M 
DO  308  I = I . N 

WFRTXI I )*WFRTX< I )+SI ZVEU#DECT  J M 
XDIFFX I  I )* WFRTXI I J -XI  I) 

IF  (  XDIFFX! I )  )  306.306.307 

306  Tin  »o, 

GO  TO  308 

307  T(  n-XCIFFXI  1  J/SIZVEL 

308  CONTINUE 
RETURN 
END 
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0! MENS JON  A NO  COMMON  STATEMENTS  ON  FIRST  PACE  OF 
COMPUTER  PRINT-OWT  ARE  TO  BE  ADOEO  HERE* 

SUBROUTINE  YIELD 
DO  77  I«1,NN 
DO  77  I  !■! *NF!BRf 

IF  (  TFEPSB< I . Itl+CEPSY  »  68*68*77 

68  IF  <  PLSETB< I , l  »- < TFEPSB < l *  1 1 I+CEPSY) I  77*77.69 

69  PLSETB (1*11) -TFEPSa (  I  ,  l  l  I ♦CEPSY 
77  CONTINUE 

DO  222  !M*NN 
R»  ABSF  <  ESEPBT ( { » ) 

IF  < R-SEPSY |  222*222(218 

218  IF  ( ESEPBT ( I ) )  219*220.221 

219  SPSTflTl  n.TSEPBTf  D+SEPSY 
GO  TO  222 

220  SPSTBTin-SPSTBTU) 

GO  TO  222 

221  SPSTBTI I»«TSEPBT< I i-SEPSY 

222  CONTINUE 

00  642  I >1 »NN 
RaABSF (ESEPBB( I ) ) 

IF  <  R-SEPSY1  642*642.638 

638  IF  (  ESEPBBlln  639.640*641 

639  SP  STBS  (11  ■  TSEPBBU)  +  SEPSY 
GO  TO  642 

640  SPSTBB( I l»SPS”BB( I) 

GO  TO  642' 

641  SPSTBB< I >«TSEPBB< 1 I-SEPSY 

642  CONTINUE 

DO  91  U2«NN 

00  91  1  I  •*  1  .NFIBRE 

IF  (  TFEPSJ ( I  *  1 1 I+CEPSYl  88*88*91 

88  IF  (  PLSETJt l *11 >-( tfepsj { I «II)+CEPSY) »  91*91*89 

89  PLSETJI I « 1 1 )«tFEPSJ( 1  *  1 ? I+CEPSY 
91  CONTINUE 

DO  262  1«2*NN 
R*» ABSF ( E5EPJT (  I)  ) 

IF  (  R-SEPSY.  262*262*298 

298  IF  (  ESEPJT ( I ) )  299*260*261 

299  SPSTjT ( I ) aSEPSY+TSEPJT (  1 ) 

GO  TO  262 

260  SPSTJTI I laSPSTjTt I ) 

GO  TO  262 

261  SPSTJT(I)«TSEPJTI I I-SEPSY 

262  CONTINUE 

DO  282  I *2 ♦ NN 
RaABSFIESEPJfcM  I ) ) 

IF  (  R-SEPSY  I  282.282,278 

278  IF  (ESEPjBdll  279.280.281 

279  SPSTJB ( 1 >»TSEPJBi I l+SEPSY 
GO  TO  282 

280  SPSTJB ( I) aSPSTJBI I) 

GO  TO  282 

281  SPSTJB! I )»TSEPJ8( I I-SEPSY 
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282  CONTINUE 

DO  78  t=l.NN 

IF  (  TEPBAHI  I  MCEPSY)  71. 

71  IF  (  EPSETUI  I  >-(  fEP0AR  I  I  ) 

72  EP'.ETU  <  i  >  =  TEPBAR  (  |  )  +CEPSY 
78  CONTINUE 

DO  92  1=2, NN 

IF  <  TEPJT{  n+CEPSY.t  20.2 

20  I F  (  EPSETJ<  I  >-<TEPJT{  I)  + 

21  EPSETJ ( I  I »TEPJT (  I ) ♦CEPSY 
92  CONTINUE 

DO  232  I =1 ,NN 
P* AbSF ( EEsEPUf I  1 ) 

IF  {  P-SKOSY)  232.232,318 

318  IF  (  EFSEPBdl  )  319.320. 

319  ESP3T8C I )*TEPBAR( I 1+SEPSY 
GO  TO  232 

320  ESPST8( I )«ESPSTBI I ) 

GO  TO  232 

321  E?,PSTB<  I  )  =Ti:P9AR<  I  1-SEPSY 
2.12  CONTINUE 

DC  272  1=2. NN 
R“ AbSF ( EESEPJ  (  I  )  ) 

IF  (  R-S£PSy)  272,272.358 

358  IF  <  EESEPJ(I)  J  399.360* 

359  ESPSTJf I )=SEPSY+T£PJT{ | ) 
GO  TO  272 

360  ESP3TJC I )=ESP5TJ(I) 

GO  TO  272 

361  ESPSTul I >»TEPJY< I I-SEPSY 
272  CONTINUE 

PE  TURN 
END 


2 .  Summary  of  FORTRAN  Notation  Used; 


The  following  is  an  alphabetic  listing  of  the  FORTRAN  notation  that 
was  used  in  the  preceding  computer  program. 


A1 

A2 

ALERB 


ALERR 

ALFA  (I) 
AQ 

AS 

AX 

AV 

B 

BARL 

BCOUNT 


-  a  pressure  read  into  the  computer  to  define  the  r ressure-t ime 
diagram;  see  Figure  11 

-  a  pressure  read  into  the  computer  to  define  the  pressure-time 
diagram;  see  Figure  11 

-  a  constant  which  when  multiplied  by  the  calculated  acceleration 
of  a-mass  gives  the  allowable  difference  between  the  assumed 
and  calculated  acceleration  of  that  mass 


-  allowable  difference  between  the  thrust  in  a  bar  or  joint 
calculated  by  r/EF0RM  and  that  calculated  by  EXTERN 

-  angle  bar  I  makes  with  horizontal 

-  absolute  value  of  the  allowable  difference  between  the  assumed 
and  calculated  acceleratior  of  a  mass 

-  total  area  of  reinforcing  .steel,  ln2/in 

-  absolute  value  of  the  difference  between  th;  calculated  and 
assumed  acceleration  in  the  x  direction  of  a  mass 

-  absolute  value  of  the  difference  between  tne  calculated  and 
assumed  acceleration  in  the  y  direction  of  a  mass 

-  Beta,  a  constant  used  in  the  Beta  Method 

-  the  original  length  of  the  bars.  Inches 

-  the  number  of  cycles  used  in  the  convergence  of  the  Beta  Method 


BLASTl(I)  -  horizontal  blast  force  acting  on  mass  I  *nd  exerting  force  to 
the  right 

BLAST2(I)  -  horizontal  blast  force  acting  on  mass  I  and  exerting  force  to 
the  left 


BLASTY(I) 

BLAST 

BMOMB(I) 

BMOMj(l) 

CEPSCR 


-  vertical  blast  force  acting  on  mass  I  and  exerting  force  downward 

-  pressure  amplitude  read  from  pressure-time  diagram  stored  in 
computer 

-  total  bending  moment  at  the  center  of  bar  I 

-  total  bending  moment  at  joint  I 

-  end  of  the  plastic  portion  of  the  concrete  stress-strain  curve; 
see  Figure  8 
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CEPSO  -  intersection  of  cr  ■  0  and  c  on  concrete  stress-strain  curve; 
see  Kigure  8 

CEPSU  -  strain  at  which  concrete  crushes;  see  Figure  8 

CEPSY  -  yield  strain  of  concrete;  see  Figure  8 

CFORCB(l)  -  thrust  in  bar  I  developed  by  concrete  h 

CFORCJ(l)  -  thrust  in  joint  I  developed  by  concrete 

CHKYl(l)  -  counter  to  determine  if  spring  1  in  the  x  direction  at  mass  I 
is  acting 

CHK Y2 ( I )  -  counter  to  determine  if  spring  2  in  the  x  direction  at  mass  I 

is  acting 

CMASS(I)  -  mass  concentrated  at  joint  I 

CMOMB(l)  -  bending  moment  at  the  center  of  bar  1  developed  by  concrete 

CMOMj(l)  -  bending  moment  at  joint  I  developed  by  concrete 

COUNTB(I)  -  the  number  of  cycles  used  in  the  correction  of  the  difference 
between  the  thrust  in  bar  I  calculated  by  DEFORM  and  that 
calculated  by  EXTERN 

CRVBAR(I)  -  the  curvature  at  the  center  of  bar  I 

CRVJT(I)  -  the  curvature  at  joint  I  •••  ,  . 

D  -  total  thickness  of  structural  element,  Inches 

DALFA(I)  -  change  in  angle  which  bar  I  makes  with  the  horizontal 

DANJT(I)  -  change  in  central  angle  at  joint  I 

DBARL(l)  -  change  in  the  length  of  bar  I 

DDXEND(l)  -  acceleration  in  the  x  direction  of  mass  I  at  the  end  of  a  time 
interval 

DDX(I)  -  acceleration  in  the  x  direction  of  mass  I  at  the  beginning  of  a 
time  interval 

ODY ( I )  -  acceleration  in  the  y  direction  of  mass  I  at  the  beginning  of  a 

time  interval 

DOYEND(I)  -  acceleration  in  the  y  direction  of  mass  I  at  the  end  of  4  time 
i nterval 

DE FX 1(1)  -  deformation  of  soil  spring  1  at  mass  I 

DEFX2(I)  -  deformation  of  soil  spring  2  at  mass  I 
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DEFV(I) 

OEFYY(l) 

DELTIM 

DEIX(I) 

DELY(I) 

{  DIFEPB(I) 

DIFEPJ(I) 

DPRIM 

r 

j  DUMMY 

DX(I) 

OXEND(l) 

DY(I) 

DYEND(l) 

i 

f 

E 

EC 

EESEPB(I) 

EESEPJ(l) 

EFORCB(l) 

EFORCJ(I) 

EFORSB(I) 

EFORSJ(I) 

EPSETB(l) 


-  deformation  of  the  vertical  soil  spring  at  mass  I 


-  yield  deformation  of  the  soil  springs  at  mass  I,  inches 

-  time  Increment 

-  displacement  of  mass  I  In  x  direction 

-  displacement  of  mass  I  In  y  direction 

-  average  strain  in  bar  I  ‘  j  correct  the  difference  between 
thrusts  calculated  by  DEFORM  and  EXTERN 

-  average  strain  in  joint  I  to  correct  the  difference  between 
thrusts  calculated  by  DEFORM  and  EXTERN 


-  distance  from  the  canter  of  the  rainforcing  ateal  to  the 
plastic  centroid 

-  a  trial  stress  computed  to  determine  stress  In  concrete 
beyond  plastic  portion  of  stress-strain  curve 

-  velocity  of  mass  I  In  x  direction 

-  velocity  of  mass  I  in  x  direction  at  the  end  of  a  time 

Interval  * 

-  velocity  of  mass  I  In  y  direction 

-  velocity  of  mass  I  In  the  y  direction  at  the  end  pf  a  tlw# 
interval 

-  dummy  variable  used  in  computer  printout 

-  concrete  modulus  of  elasticity,  pel 

-  average  "effective"  strain  In  bar  I;  average  strain  In  bar  I 
minus  the  plastic  set 

-  average  "effective"  strain  In  joint  I;  average  strain  in 
joint  I  minus  the  plastic  set 

-  concrete  force  In  bar  I  computed  using  average  strain  for  1 

that  bar 

-  concrete  force  in  joint  1  computed  using  average  strain  for 
that  joint 

-  steel  force  in  bar  I  computed  using  average  strain  for  that  bar 

# 

-  steel  force  in  joint  I  computed  using  average  strain  for  the 
joint 

'M 

-  plastic  set  in  bar  I  arising  from  the  use  of  average  strain  for 
that  bar 
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EPSETJ(I)  -  plastic  set  in  joint  I  arising  from  the  use  of  average  strain  for 
that  joint  ’  '  ■  ; 

EPSJT(l)  -  average  strain  at  joint  I 

EPSNEG  -  the  number  of  concrete  fibers  which  have  a  negative  strain  at 
time  t 

'f  '  1 

ES  -  steel  modulus  of  elasticity,  psi 

ESEPBB(l)  -  "effective"  strain  in  the  bottom  reinforcement  of  bar  I; 
bottom  steel  strain  minus  plastic  set 

!  ',X 

ESEPBT(l)  -  "effective"  strain  in  the  top  reinforcement  of  bar  I;  top 
steel  strain  minus  plastic  set 

ESLPJB(l)  -  "effective"  strain  in  the  bottom  reinforcement  of  joint  J; 
bottom  steel  strain  minus  plastic  set 

ESEPJT(I)  -  "effective"  strain  in  the  top  reinforcement  of  joint  I;  top 
steel  strain  minus  plastic  set 

ESPSTB(l)  -  plastic  set  in  steel  reinforcement  in  bar  I  arising  from  use 
of  average  strain  of  that  bar 

ESPSTJ(l)  -  plastic  set  in  steel  reinforcement  in  joint  I  arising  from 
use  of  average  strain  of  that  joint 

EYEPSC  -  compressive  yield  strain  of  concrete  using  average  bar  and 
joint  strains;  CEPSY  minus  plastic  set 

FH1(I)  -  horizontal  component  of  thrust  in  bar  I  at  end  nearest  joint  1 

FH2 ( I )  -  horizontal  component  of  thrust  in  bar  I  at  end  nearest  joint 

I  +1 

PH(l)  ”  dual  notation;  horizontal  component  of  thrust  in  bar  I;  total 
horizontal  force  on  joint  I 

1. 1  !.•!•  i  | 

FIBRE  -  the  number  of  concrete  fibers  used  to  represent  tne  cross- 

section;  limited  to  a  maximum  of  20  J! 

FORCB(I)  -  total  force  in  bar  I 

FORCj(l)  -  total  force  in  Joint  I 

FORCY(I)  -  yield  forces  in  soil  springs  at  mass  I,  pel 

FORDFB(l)  -  difference  In  *he  thrust  in  bar  I  calculated  by  DEFORM  and 
that  calculated  by  EXTERN 

FORDFj(l)  -  difference  in  the  thrust  in  jolt;  I  calculated  by  DEFORM  and 
that  calculated  by  EXTERN 
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FSIGB (1 , II)  -  concrete  stress  in  fiber  II  of  bar  I 

FSHJd.II)  -  concrete  stress  in  fiber  II  of  Joint  I 

FVl(l)  -  vertical  component  of  thrust  in  bar  I  at  end  nearest  joint  I 

FV2 ( I)  -  vertical  component  of  thrust  In  bar  I  at  end  nearest  joint  I  +  I 

FV(I)  -  dual  notation;  vertical  component  of  thrust  in  bar  1;  total 

vertical  force  on  joint  I 

FX1(I)  -  horizontal  force  on  mass  I  from  soil  spring  1 

FX2(I)  -  horizontal  force  on  mass  I  from  soil  spring  2 

FY(I)  -  vertical  force  on  mass  I  from  vertical  soil  spring 

GAMMAS  -  weight  of  a  cubic  inch  of  soil,  pci 

HGT  -  a  constant  which  when  multiplied  by  a  bar  length  gives  radial 

distance  to  which  soil  is  assumed  to  respond  with  structure 

JC(I)  -  counter  used  In  the  determination  of  structural  collapse 

jY ( I)  -  counter  used  in  the  determination  of  structural  collapse 

KOUNT  -  a  counter  to  control  printout  of  the  computer  program 

Lc (I)  -  counter  used  in  the  determination  of  structural  cdllapse 

LY(I)  -  counter  used  in  the  determination  of  structural  collapse 

N  -  total  number  of  joints  in  structure;  limited  to  a  maximum  of  11 

NFIBAE  -  number  of  fibers  used  to  represent  the  concrete  cross  section; 

limited  to  a  maximum  of  20 

NN  -  total  number  of  joints  minus  one;  number  of  bars  in  structure 

NPRINT  -  a  constant  read  into  computer  to  control  the  printout,  usually  *  2 

PALFA(I)  -  angle  bar  I  makes  with  the  horizontal  at  the  beginning  of  a 
time  Interval,  degrees 

PDDX(l)  -  the  assumed  acceleration  of  mass  I  In  x  direction  for  any 
trial  iteration  in  the  Beta  Method 

PDDY(I)  -  the  assumed  acceleration  of  mass  I  in  y  direction  for  any 

trial  iteration  in  the  Beta  Method 

PHEE  -  angle  wave  front  makes  with  the  horizontal 

PLSETB(l ,11)  -  plastic  set  in  concrete  fiber  II  of  bar  I 

PLSETJ(I.II)  -  plastic  set  in  concrete  fiber  II  of  joint  I 
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px(!)  -  x  coordinate  of  mass  I  at  the  beginning  of  a  time  Interval 

PY(*)  “  y  coordinate  of  mass  I  at  the  beginning  of  a  time  interval 

Q  -  constant  used  in  the  location  of  the  wave  front 

Oft  “  absolute  difference  between  thrusts  calculated  by  DEFORM  and 

those  calculated  by  EXTERN 

R  -  a  constant  used  to  determine  the  time  history  of  strains  In 

concrete  and  steel 

RATIO  “  ratio,  of  horizontal  to  vertical  blast  pressures  on  masses; 

read  into  machine 

SCRATO  -  ratio  of  ES  to  EC  '  ! 

SEPSY  -  yield  strain  in  steel  reinforcement 

SEPYPC(I)  -  compressive  yield  strain  in  steel;  plastic  set  minus  SEP^Y 

SEPYPT(I)  -  tensile  yield  strain  in  steel;  plastic  -  et  plus  SEPSY 

SFORBB(l)  -  force  in  the  bottom  reinforcing  steel  of  bar  I 

SFORBT(l)  -  force  In  the  top  reinforcing  steel  of  bar  I 

SFORJB(I)  -  force  In  the  bottom  reinforcing  steel  of  joint  I 

SFQRJT(I)  -  force  In  the  top  reinforcing  steel  of  joint  I 

SIGY  -  dummy  variable  used  In  computer  printout 

SIZVEL  -  velocity  at  which  wave  traverses  the  structure,  in/sac 

SPSTBB(l)  -  plastic  set  in  the  bottom  reinforcing  steel  of  bar  I 

SPSTBT(l)  -  plastic  set  in  the  top  reinforcing  steel  of  bar  I 

SPSTJB(I)  -  plastic  set  in  the  bottom  reinforcing  steel  pf  Joint  I 

SPSTJT(l)  -  plastic  set  in  the  top  reinforcing  steel  of  Joint  I 

T(I)  -  total  time  the  wave  has  been  in  contact  with  mass  I;  used 

to  determine  pressures  from  a  pressure-time  diagram 

TAR2  -  time  of  arrival  of  A2;  see  Figure  11 

TDELX(I)  -  total  displacement  of  mass  I  in  x  direction 

TDELY(I)  -  total  displacement  of  mass  I  in  y  direction 

TEPBAR(l)  -  total  average  strain  in  bar  I 

TFEPSBd.H)  -  total  strain  in  concrete  fiber  II  of  bar  I 
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TFEPSJd.II)  -  total  strain  trt  concrete  fiber  II  of  joint  I 


TFINAL 


TSEPBB(I) 

TSEPBT(I) 


time  at  which  computer  Is  shut  off 

total  elaspsed  time  since  wave  has  contacted  structure 

time  required  for  blast  to  rise  to  pressure  Of  Al;  see 
Figure  11 

time  at  which  pressure  attains  an  amplitude  of  A2 ;  see 
Figure  11 

total  strain  in  the  bottom  reinforcing  steel  of  bar  I 
total  strain  In  the  top  reinforcing  steel  of  bar  I 


TSEPJB(l)  -  total  strain  In  the  bottom  reinforcing  steel  of  Joint  I 

TSEPJT(l)  -  total  strain  in  the  top  reinforcing  steel  of  joint  I 

VB ( I)  -  shear  at  the  ends  of  bar  I 

VH(I)  -  horizontal  component  of  shear  In  bar  I 

VH 1(1)  -  horizontal  component  of  shear  In  bar  I  at  end  nearest  Joint  I 

VH2(I)  -  horizontal  component  of  shear  in  bar  I  at  end  nearest  joint 

I  +1 

VV  -  dummy  variable  to  determine  the  presence  of  horizontal  iprfnjjjs 

1  and  2 

VV(l)  -  vertical  component  of  shear  In  bar  I 

VVl(l)  -  vertical  component  of  shear  In  bar  I  at  end  nearest  Joint  I 

VV2(I)  -  vertical  component  of  shear  In  bar  I  at  end  nearest  Joint 

I  +  1 

WFRTX(l)  -  location  of  wave  front  with  respect  to  joint  I  of  structure 

X(I)  -  x  coordinate  of  mass  I  at  any  time 

XDIFFX(I)  -  distance  between  wave  front  and  mass  I  of  structure 

XORIG(I)  -  original  x  coordinate  of  mass  I 

Y(I)  -  y  coordinate  of  mass  I  at  any  time 

YC  -  a  constant  used  in  the  determination  of  structural  collapse 

YCOUNT  -  a  constant  specifying  failure  criteria  for  structure 

YEPSC  -  a  compressive  yieH  strain  of  concrete;  CEP$Y  minus  plastic  fft 
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YORIG(I)  -  original  y  coordlnatt  of  mass  I  ..  *.  ! *  i'l1''’ 'i  j 

-  dummy  variable  to  determine  the  presence  of  horltontdl  iprinqs 

1  and  2  1  :,K  ffl&jj 

;  1  !  . . 11  " 1  ''  1  »  1  ,  *■; { ■! vl:i!iiij(  % ii|„iu  i  ■ 

i*  variable  to  determine  the  presence  Of  her liontal  spring*  . 

1  #nd  2  '  1  ■  ::::|p5ip:! 

,,  ,  *  1  •:*,!;  1 , h'iiiji;  'i  '5i  I j't'l 1 

-  dummy  variable  to  determine  the  presence  of  horizontal  springs 
1  and  2 


wr  noriionmi  springs 

•  !M'!  ’ 

- 1  »ii  r/y»|i 


,1:,  'I1  ,  ,«'■ 


II''  hi 


jiW!  ;|f« ,  .  IN' 

ffl  in  • 1  I U 


'•  Pm1 

■  i.  v  % 


"  'rt 


: 

:i!  'I  ’ 

•  I  !l 


m\\ 

: 'r!i:|'  I!'!,  phi 

'•'I  "  !  fVjI'i, 


1)1  1  Vwi: 


-i  4  i  <!•  m 


,;l  if  :ji  ’'i  '!i| 


'i  M,  X  ,h  i  il  i:  1  i|r  |  j 


'ill11''  'iiii  i'i", 


" . . . . . il  t  I'n  'inriii’if  "nmrmnit  . . 


! II  r^:. it;  il  i;  JI I1U  LiLiiMM!  IMM  j  jjl!  H  ^1 1  [lilp  imN 


inWfW'M 


'  '1  •s 


3-  iMnnarv  of  Inou £  &»*«• 


■t  I'll'i  (1  I') 


data  Iiput°[i°^he9cOTDlterTJo«ra!y  “il*  an?  foriWt'  of  the  b«ic 

_  p  ,  ,  ,  tne  computer  program.  The  data  are  listed  bv  their  FOR tram 

"***•  wh,ch  *»  Identlfl.d  by  reference  to  Sub-t.ctl^  2  " 


Card  No. 

1 

2 

3 

4  +  a* 

5  +  2a 

6  +  2a 

7  +  2a 

8  +  2a 

9  +  2a 

10  +  2a 

11  +  2a 

12  +  2a  +  b** 

13  +  2a  +  b 

14  +  2a  +  b 

15  +  2a  +  b 

16  +  2a  +  b 


FORMAT 

13.E12.4 

5E12.4 

E12.4 

5612.4 

5E12.4 

E12.4 

2E12.4 

I3.E12.4 

3E12.4 

5E12.4 

2E12.4 

5E12.4 

2E12.4 

E12.4 

2E12.4 

13 


mm  am 

N,  BARE 

A1 ,A2,TR1 ,TAR2,  TR2 
RATIO 

FORCY(I)  I  -  2.  N-l 
DEFYY(I)  I  .  2,  N-1 
S1ZVEL 

DELTIM,  TFINAL 

nfibre,alerr 

D.DPRIM.A5 

ECtCEPSY,CEPSCR.:r,PSU,CEPSO 

E$,$EP$Y 

PALFA(I)  I  -  1,  N-l 

B.ALERB 

PHEE 

GAMMAS, HOT 
NPRINT 


*  If  N-l  >  6  a  second  card  will  be  required. 

**  If  N-l  >  5  a  second  card  will  be  required. 


